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ABSTRACT 


Two recent computer programs, WRCOMP and ENGINE, by A. 
Mathur, were used to examine the performance to be gained by 
incorporating a wave rotor component in a turbofan engine with 
mixed exhausts. The programs were transferred to a VAK-2009 
computer, extended, and test cases reported by A. Mathur were 
successfully reproduced. A comparison was made between 
ENGINE, in which real gas effects are accounted for, and ONX 
(by J. Mattingly) in which constant specific heats are used. 
The inclusion of real gas effects proved to have a significant 
impact on the predicted performance. An extension of Mathur's 
results was made by varying the overall pressure ratio in the 
wave-turbofan engine. Further cycle studies and experiments 
to measure wave rotor component performance are 


recommended. 
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I. INTRODUCTION 


Interest by DARPA (Defense Advance Research Projects 
Agency) and NAVAIRSYSCOM (Naval Air Systems Command) in 
obtaining increased range for unmanned air vehicles (UAV's) 
prompted the present study of wave rotor applications in 
engines. The wave rotor offers potential to achieve lowes 
specific fuel consumption (SFC) and higher specific thrust by 
allowing increased cycle temperatures without increasing 
turbine inlet temperature. The goal of this study was to 
examine the potential benefits of incorporating a wave rotor 
in a turbofan engine with mixed exhausts. The results showed 
that significant benefits could be obtained in both SFC and 
specific thrust denending upon the selection of other cycle 
parameters. 

Chapter II describes the wave rotor concept, presents a 
general background on trends in gas turbine technology, gives 
a description of a wave rotor engine component and introduces 
the two computer programs (ENGINE and ONX) which were used in 
the study. Chapter III gives a detailed description of the 
ENGINE code and presents a comparison of results obtained with 
ENGINE and ONX. Chapter IV gives a description of the WRCOMP 
program, which analyses the unsteady wave rotor flow, and its 
interface with the ENGINE program. Chapter V gives results 


obtained using the codes for the performance of a turbofan 








with mixed exhausts, with and without a wave rotor component. 
Conclusions and recommendations for further work are given in 
Chapter VI. Details are given in the appendices. Appendices 
A and B contain procedural guides, sample results and program 
listings for ENGINE and WRCOMP, respectively, while Appendix 
C describes file transfer and graphics procedures which were 


used in the course of the study. 





II. BACKGROUND 


A. POTENTIAL FOR WAVE ROTOR APPLICATIONS IN ENGINES 

A wave rotor is a partial admission device which causes 
one gas to compress another by wave propagation [Ref. l:p. 1]. 
A series of straight compression tubes are aligned on a rotor 
encased in adrum. The rotor rotates on a fixed axis. Flow 
in and out of the wave rotor is steady while flow inside the 
wave rotor tubes is unsteady. Figure 1 shows an example of a 
wave rotur component. The single wave rotor component serves 
the same function as the combination of a compressor and a 


turbine mounted on the same shaft. 


Compressed alr 





Figure 1. Wave Rotor Component [Ref. l:p. 294] 


Figure 2 shows a turbofan engine configured with an 


axially mounted wave rotor component. 











Figure 2. Wave Rotor Turbofan Engine [Ref. l:p. 293] 


Development of wave rotor devices for engines began in the 
early 1950's. By the early 1960's, gains were achieved in 
performance through improved seal design and improved timing 
of the waves arriving at inlet and outlet ports ("tuning") 
{Ref. l:p 174}. Few developments were attempted after the 
1960's, except by Brown-Boveri in Switzerland. The Brown- 
Boveri "“Comprex" is presently being produced as a commercial 
supercharger and is available in trucks, heavy equipment and 
Mazda diesel automobiles. Reference 1 provides a comprehen- 
sive review of wave rotor technology. 

Wave rotor components have promising applications in small 
engines in the 600 to 1000 lb. thrust range. Such engines 
have applications in cruise missiles, remotely piloted 


vehicles (RPV's), helicopters and small thrust turbofans. The 


more conventional epproach, however, has been to pursue 


increased turbine inlet temperatures in conventional gas 











turbine engines to achieve the desired improved performance. 
However, in the process, the requirement for turbine cooling 
also increased. The wave rotor offers the potential to 
increase performance without the attendant necessity for 
cooling. 

Studies of engine performance indicate that a wave rotor 
component must be included in the initial cycle optimization 
rather than be added to an existing engine design. Wave rotor 
components require, optimally, different compressor and fan 
pressure ratios. The inclusion of a wave rotor reduces the 
number of stages required for compression to a constant cycle 
pressure ratio and allows higher maximum temperatures in the 
cycle for the same turbine inlet temperature. Performance 
predictions for gas turbine engines incorporating wave rotor 
components show decreased specific fuel consumption and 
increased specific thrust for any fixed value of turbine inlet 
temperatures. 

Figure 3 shows the gains that are thought to be achievable 
in the performance of small engines [Ref. l:p. 292]. In Figure 
3, “near term improvements" include higher rpm (for a given 
engine size) and higher turbine temperatures. "New engines" 
will include composites, ceramics, or new materials. Some 
current "revolutionary engines" include the compound cycle, 
the eccentric, the recuperated and the wave rotor engine [Ref. 
l:p. 292}. The conventional, eccentric and the recuperated 


engines are fundamentally limited by the maximum allowable 











TSFC (IbAb-ty) normalized to currant engines 





Figure 3. Trends in Engine Development [Ref. l:p. 292] 


turbine inlet temperature (TIT) [Ref. 2:p. 60]. The wave 
rotor engine is the only engine that promises significantly 
increased performance without requiring an increase in TIT 


(Ref. l:p. 293}. 


B. ENGINE PERFORMANCE CODES 

Two engine performance codes were available for the 
present study, ENGINE and ONX. There were some differences 
between the two programs. ENGINE calculated real gas effects 
based on fuel-air ratio and temperature, while ONX used 


constant specific heats at the compressor outlet and 








downstream of the turbine inlet station. In both programs, 
the engine performance is obtained by equating power balances 
across the turbine and compressor sections, and describing an 
expression for the enthalpy rise a ~oss the burner section. 

Program ENGINE was developed by A. Mathur [Refs. 3; 4], 
and uses notation developed by Vavra [Ref. 5], and Vanco 
{Ref. 6]. The program calculates the performance of a 
turbofan engine with mixed exhausts with and without the wave 
rotor component. 

Program ONX was developed by J.D. Mattingly [Ref. 7], and 
it uses methods and notation developed by G. C. Oates [Ref. 
8]. The program calculates the performance for a variety of 
engine configurations including a turbofan with mixed 


exhausts. 


C. WAVE ROTOR INTERNAL FLOW DESCRIPTION 

In 1984, A. Mathur developed a FORTRAN program (WRCOMP) to 
analyze the unsteady processes in wave rotors [Ref. 9]. A 
detailed description of the analysis and code are provided in 
Reference 9. WRCOMP uses a 1-D, random choice method (RCM) to 
solve the hyperbolic, nonlinear conservative system of 
equations by solving localized Riemann problems. Shocks have 
sharp resolution. Discontinuities are modelled as jumps. The 
code is first order accurate in time and it is unconditionally 


stable. Further explanations and details on the Random Choice 





Method are given by Glimm [Ref. 10], Chorin [Ref. 11], and Sod 
[Ref. 12]. 

WRCOMP calculates the unsteady process inside the wave 
rotor, outputs the inlet and outlet port opening and closing 
times required for correct design, and computes the perfor- 
mance of the wave rotor given the total pressure ratio (TPR), 
the static pressure ratio (SPR), the exit Mach number (AMEXIT) 
and the (initial) reference density inside the wave rotor. 
The performance calculations depend on a favorable comparison 
for the mass flow rate into and out of the wave rotor 
component. The outputs from WRCOMP are used as inputs to 
engine performance calculations performed using the ENGINE 


code. 








III. ENGINE PERFORMANCE CODE 


A. DESCRIPTION OF THE CODE 

A complete listing, procedural guide and sample results 
for ENGINE are provided in Appendix A. The code incorporates 
modular subroutines to calculate each section of the engine 
cycle. A complete derivation of the code can be found in 
Reference 4. 

The following changes were made to improve the utility of 
the code: 


(1) The original data statement format was modified for a 
user input interface; 


(2) A pop-up menu depicting the configurations was added 
(see Appendix C for details); 


(3) An altitude table was incorporated to automate ambient 
conditions using Reference 13; 


(4) A graphics plotting routine was added using GRAFkit 
software, Reference 14; 


(5) A loop was constructed to increment for: 
a) bypass ratio, 
b) fan pressure ratio, 
c) compressor pressure ratio, 
ad) LPT exit Mach. 
A complete listing of the input parameters and a sample output 
result appears in Appendix A. 
ENGINE (and ONX) computes engine performance in terms of 


uninstalled specific thrust and specific fuel consumption. 




















The uninstalled specific thrust is defined as the thrust per 
unit mass flux through the engine and determines the encine 
size. The specific fuel consumption is the mass flow rate of 
fuel per unit thrust and determines the range. A power 
balance is written for each turbine and compressor spool. An 
enthalpy rise is written for the flow through the combustor. 
ENGINE uses loss coefficients and efficiencies to account for 
component performance. The program also uses real gas effects 
based on fuel air ratios and temperature. In contrast, ONX 
uses pressure ratios and polytropic efficiencies to account 
for component performance and assumes constant specific heais 
across the compressor (Cp,) and downstream of the combustor 
(Cp.). A necessary step, prior to the use of ENGINE, or 
comparison between ENGINE and ONX, was to reproduce the 
results reported by A. Mathur in Reference 4. These results 
were successfully duplicated using the inputs given in 
Reference 4 and Appendix A. 

To achieve an accurate comparison between ENGINE and ONX, 
the polytropic efficiencies required in ONX were computed from 
the outputs of the ENGINE code using the component efficien- 
cies, the actual and isentropic pressures and temperatures, 


and gamma. A sample calculation is provided in Appendix A. 


B. COMPARISON WITH ONX CODE 
Since engine cycle calculations are greatly affected by 


the inputs, especially the component efficiencies, an accurate 








comparison between ENGINE and ONX could only be performed if 
the flight parameters, the design choices, the design 
limitations and the component efficiencies were the same. 
Flight parameters include the flight Mach number and the 
altitude. Design choices include the maximum temperature in 
the burner, the lower fuel heating value, the maximum turbine 
inlet temperature, the low pressure turbine exit Mach number, 
the compressor pressure ratio, the fan pressure ratio, the 
bypass ratio, etc. Fans, compressors, burners, turbines and 
nozzles are characterized by component efficiencies, which are 
the coupline:.ts to the loss coefficients. The simplest test 
case for comparison of ENGINE and ONX was the non-afterburning 
turbojet. ENGINE calculated the turbojet case by setting the 
bypass ratio equal to zero. 
1. TURBOJET 

With all input and component performance parameters 
matched, the results of ENGINE were compared to those of ONX 
for the case of the turbojet. The results are shown in Table 
I. For the ONX (1) case, the Cp, and Cp, were selected to be 
in the range that Mattingly referred to as "typical" values 
(Ref. 7:pp. 116, 438]. The ONX (2) values for Cp, and Cp, were 
the result of numerous iterations using ONX to try to achieve 
agreement with ENGINE. The ONX (3) values for Cp, and Cp, were 


obtained by computing the average values for Cp, and Cp, from 


11 





TABLE I 
TURBOJET PERFORMANCE COMPARISON 


M, = 0.79, h=35000 ft. 


ENGINE: F/m 
SFC 


41.88 lbf/(lbm/sec) 
0.899 (lbm fuel/hr)/lbf thrust 


ONX ONX ONX 


(1) (2) (3) 


Typical Values req'd Average 


Specific for from ENGINE 
agreement results 


Difference 
F/m 


Difference 
SFC 





the ENGINE program output data. The differences were seen to 


be surprisingly large 


A comparison table was constructed of the pressure and 


temperature into each component to trace where’ the 


calculations began to diverge. The divergence was found to 


occur in the enthalpy rise calculation across the combustor. 


For ONX, the enthalpy rise gives the fuel-air ratio (f) as 


12 


Cp,xT,4~CD-xT,,; 


ee 1 
[Cp,xT,.+h,xn »] -Cp_xT yz, om 


and the test case in Table I gave f = 0.0152. For ENGINE, the 


fuel-air ratio equation was: 


va A(T,,) -A(T,3;) +£;x (B(t,,) -B(t,3)] (2) 
(Cp xT, .+hyxn p] -B(T,4) 


where, 
Cc. C. ‘as c 
A(T) =Cyx Toe kT? 4 xT? = ex Tee xl? (3) 
and 
Dy: D D D 
B(T)=D xT+—2xT?-—3x7T?+—4 xT*- 3 x75 4 
(T)=D_x ax ; . (4) 


With f, = 0.0 (f prior to the burner),the test case designated 
ONX (1) in 7” le I gave f = 0.0103. Thus ONX calculated a 
32.24% higher value for the fuel air ratio than did ENGINE. 
The question then was, can agreement be obtained between ONX 
and ENGINE using any reasonable values of the Cp's in ONX. 
The result of an attempt to achieve such agreement is shown as 
ONX (2) in Table I. Using Cp, = 0.235 and Cp, = 0.251, the 
disagreement in specific thrust was reduced to 4.4% and in SFC 
to 6.1%. It is noted here that the Cp value for air at the 


temperature at the exit of the combustor (1860°R) in the 


13 





present example is Cp = 0.271, considerably higher than the 
value required to approach agreement. 

One further attempt was made to compare the code 
predictions for a turbojet. Values of Cp were calculated at 
the inlet and outlet stations of the compressors and turbines 
using the enthalpies calculated at the stations by ENGINE in 
the expression, h = CpT. An average of the in and out values 
was then used as an input to ONX. The results are shown as 
ONX (3) in Table I. The strong sensitivity of the performance 
to the values used for the specific heats is evident in the 
comparison of the results of ONX (3) with ONX (2). 

The selection of Cp values, therefore, had a 
significant affect on the performance calculations for 
specific thrust and specific fuel consumption. In order to 
examine the sensitivity to specific heats, ONX was run with 
varying Cp, and varying Cp,. Figures 4 and 5 show the affect 
that Cp, and Cp, had on the specific thrust and SFC given by 
the ONX code. The sensitivity to the selection of the input 
values for ONX is clearly evident in these figures. 

A complete deviation of the real gas equations 
programmed in ENGINE can be found in Reference 3. Further 
details can be found in References 5 and 6. In short, the 
stoichiometric combustion equation for a general hydrocarbon 


fuel (CH,), is expressed as a polynomial expansion. The 


constants in the expansion C, to C, and D, to D, in Equations 
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Figure 4. Effect of Cp, on Specific Thrust and SFC 


Cyie G28 yee 43850, Cpn Ss 0.8434) 


(3) and (4) are deiived from the constituents of the combus- 
tion gases. The final form of the enthalpy rise for the 
combustor is given as Equation (2). Thus Cp is effectively, 
and properly, a function of fuel-air ratio and temperature. 
Since the ENGINE code calculated the enthalpy using the fuel- 


air ratio and temperature rather than maintaining constant 
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Figure 5. Effect of Cp, on Specific Thrust and SFC 


(¥.= 1.400, ¥, = 1.350, Cp, = 0.2508) 


values for specific heats, the ENGINE code results were now 


considered to be more reliable than the ONX results. 


2. TURBOFAN 


After successfully reproducing A. Mathur's results 
[Ref. 4:p.7] for a turbofan with mixed exhausts, a comparison 


between ENGINE and ONX was carried out for the same engine. 
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input into ONX. The results are presented in Table II. The 
percent differences in the specific thrust and SFC, predicted 


by ONX in comparison to ENGINE were 11% aid 23% respectively. 


[Ref. 


The average values of Cp's given by the ENGINE results were 


4:p. 7} 


TABLE ITI 


TURBOFAN PERFORMANCE COMPARISON 


M.= 0.79, h=35000 ft. 


ENGINE [Ref. 4:P.7]: 
F/m 30.185 lbf/(lbm/sec) 
SFC 0.832 (lbm fuel/hr)/1lbf thrust 


ENGINE ONX 





Cp,=0.2434 


Cp,=0.2508 


2° 


% 
Difference 
F/m 
% 
Difference 
SFC 


The bypass ratio was varied keeping other parameters 
fixed. Using ENGINE, a very narrow range of bypass ratio was 
found for which the static pressures of the bypass stream and 
the core flow could be equal at the constant pressure splitter 
plate aft of the low pressure turbine. The results in Figu.-e 
6 show that common trends were obtained with ENGINE and ONX. 
As bypass ratio was increased, both SFC and specific 
thrust decreased. However, the effect on specific thrust was 
the largest. In fact, only a 5% decrease in SFC could be 


obtained at the cost of a 20-30% decrease in specific thrust. 
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Figure 6. Engine Performance Trends for ONX and 
ENGINE Varying Bypass Ratio 


(M, = 0.79, h = 35000 ft, nm, = 2.3, n, = 10.87, TIT = 1860 °R) 
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To complete the comparison of ENGINE and ONX and to 
examine the predictions of the ENGINE program, the following 
parameters were varied and the results obtained for specific 
thrust and SFC were plotted in the indicated figures: 1) exit 
Mach number from the LPT (Figure 7); 2) maximum temperature in 
the combustor (Figure 8); 3) fan pressure ratio (Figure 9); 


and 4) compressor pressure ratio (Figure 10). 
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Figure 7. Effect of LPT Exit Mach Number (AMACH8) 
on Specific Thrust and SFC 


(M, = 0.79, h = 35,000 ft, mn, = 2.3, mn, = 10.87, TIT = 1860°R) 








As shown in Figure 7, the LPT exit Mach number had 
little effect on the specific thrust or on the SFC. 

In Figure 8, results are shown for varying turbine 
inlet temperature obtained using ENGINE with bypass ratio 
allowed to vary and with ONX with bypass ratio fixed. In the 
cases shown in Figure 8, the bypass ratio of ENGINE was 
allowed to vary such that the static pressure at the splitter 
plate, where the bypass air is mixed with core air, was 
matched. The mixed exhaust forced higher values of bypass 
ratio as the TIT was increased. The higher bypass ratio 
allowed the specific thrust to remain fairly constant. Since 
significantly more air was bypassed, less fuel was required to 
produce the same specific thrust as the TIT was increased. 
Since the bypass ratio increased, the SFC decreased as the TIT 
increased. 

The results from ONX plotted on Figure 8 show the 
effect of increasing TIT while keeping bypass ratio constant. 
Fixing the bypass ratio resulted in increased specific thrust 


but at the price of increased SFC. 
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Figure 8. Effect of Varying TIT on Specific Thrust and SFC 


(M, = 0.79, h = 35,000 ft, ny = 2.3, nm, = 10.87, TIT = 1860°R) 


The range of fan pressure ratio was taken to be 
between 1.0 and 4.0. Figure 9 shows the performance 
calculated as the fan pressure ratio was varied between 1.0 
and 3.4. ENGINE required that the bypass ratio vary to 
satisfy the equal-static-pressure condition at the splitter 
plate. A very narrow band of fan pressure ratios was obtained 
if the bypass ratio was also specified. ONX produced a 


broader range of values for a specified bypass ratio. 
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Figure 9. Effect of Varying PRLC on Specific Thrust and SFC 


(M, = 0.79, h = 35,000 ft, n, = 10.87, TIT = 1860°R) 


The result of varying compressor pressure ratio is 
shown in Figure 10. The bypass ratio was fixed at B = 0.5. 
ENGINE predicted a minimum SFC for compressor pressure ratio 
range between 11.5 and 12.5 (whereas the minimum predicted by 
ONX was greater that 20). It is clear from the curves in 
Figure 10 that a pressure ratio less than that for optimum SFC 
would be chosen to obtain a significant improvement in 


specific thrust. 
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IV. WAVE ROTOR ANALYSIS 


A. SUMMARY DESCRIPTION 

The wave rotor component and the combustion chamber can be 
considered to be combined to form a two port component called 
a gas generator. Other models of the wave rotor component are 
possible, but the gas generator model is the most straight- 
forward. The equations and code for the wave rotor flow 
modeling used in the WRCOMP program are described in 
Reference 9. The equations for the wave rotor component as a 
gas generator included in the ENGINE code are derived in 
Reference 4. Figure 11 shows an unwrapped wave rotor and a 
simplified wave diagram [Ref. l:p. 38]. Reference 1 describes 
the gas generator process as a filling, emptying, filling and 
emptying process. 

The process can equally be described as two scavenging 
processes separated by periods of stationary flow within the 
rotor. There is a low pressure scavenge at station (4) anda 
high pressure scavenge at station (2). The combustion exhaust 
gas initially compresses the air inside the wave rotor and 
then is scavenged to (4) by the incoming air at the low 
pressure port (1). When the high pressure exhaust port opens, 


the high pressure gas inside the wave rotor exits to the HPT 


and effectiv2ly scavenges the wave rotor [Ref. 15:p. 63]. 
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Figure 11. Unwrapped Wave Rotor Component [Ref. l:p. 38] 


In the ENGINE program, the gas generator is modelled as an 
internal combustion process in which a gain in total pressure 
and a rise in the total temperature is achieved. The 
performance of the gas generator is expressed as a function of 
the static pressure ratio, the total pressure ratio and the 


total temperature ratio across the inlet ports [Ref. 3:p. 14]. 


B. INTERFACE WITH ENGINE PERFORMANCE CODE 

WRCOMP can be run alone or in conjunction with the ENGINE 
code. Since the VAX system operates in windows, both programs 
can be run simultaneously. Both procedures for running WRCOMP 
are outlined in Appendix B. The original code used data 


statements in SI units. It was modified to accept input 


pressures in psi, temperatures in °R, density in lbm per ft’, 








gamma, fuel air ratio (f), and wave rotor exit Mach number. 
In the coding, the analysis is nondimensionalized to the input 
parameters. An output file, WRPERF1.0UT, was incorporated to 
store the input parameters and to convert the English units to 
SI units. Routines were added to compute the procedures 
outlined in Reference 4 and to send the output to WRPERF1.OUT. 
The DISSPLA plotting routines were disabled. In practice, 
WRCOMP is run until the mass flow into the wave rotor is equal 
to the mass flow out. Adjustments are made progressively to 
RREF, the reference density inside the wave rotor, until the 
mass flows are equal. For 1681 iterations, a CFL number = 
0.6, a grid cell width of 0.01, and performing 20 Riemann 
iterations, the total run time was about 14 minutes on the 
VAX. An example output is presented in Appendix B. 

From experience with WRCOMP and ENGINE, it is easier to 
keep the overall pressure ratio constant for a series of test 
cases. By keeping the overall pressure ratio constant, the 
outlet pressure from the HPC stays the same. If the outlet 
pressure from the HPC is the same and the performance 
parameters for the wave rotor are kept the same, then the 
value for the reference density which satisfies the mass flow 
solution for WRCOMP also stays the same. 

In practice, it is better to vary the overall pressure 
ratio and maintain the same turbine inlet temperature. This 
allows optimizing the fan, the compressor and the wave rotor 


based on realistic cycle constraints. 
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Vv. LOW BYPASS TURBOFAN CYCLE CALCULATIONS 


A. GENERAL CONSIDERATIONS 

The present work has been concentrated on the low bypass 
turbofan engine, and what performance improvements can be 
obtained by incorporating a wave rotor component. The initial 
calculations have been made without reference to a specific 
vehicle and mission. Considerations which underlie the 
calculations and which limit their scope are discussed in the 
following paragraphs. 

The characteristics of low bypass turbofans provide a 
reasonable match to the performance requirements for cruise 
missiles and RPV's. However, the design limitations for the 
components of a cruise missile may vary significantly from 
those of an RPV, even though the two engines may have very 
similar thrust requirements. A cruise missile may require 
compact packaging and "low cost design" as it is considered a 
"throw-away" or one-time-use engine, thereby relaxing the TIT 
and cooling requirements. An RPV may require a lower TIT to 
prolong its service life. 

Some limitations are inherent in the bypass engine type. 
A turbofan with mixed exhausts may be bypass ratio limited. 
The bypass air is mixed with the core air at the splitter 
plate, where the static pressures must match. If the Mach 


number of the bypass air is too low, the losses are too great. 
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Reasons for mixing the exhaust streams are to reduce engine 
noise and to reduce the IR signature by reducing the 
temperature of the core air by mixing with the cold bypass 
air. There is an improvement in performance ideally if mixing 
occurs internally when there is a large temperature difference 
between the two streams. [Ref. 8:pp. 165, 172] 

The present representation of the wave rotor component is 
somewhat limited. The modeling associated with instantaneous 
port opening and closing may require further investigation. 
Slow port openings may provide more efficient wave rotor 
operation at length to width ratios below three and stagnation 
pressure losses are predicted to increase as the length to 
width ratio increases from three to six [Ref. l:p. 236]. 

In performing cycle calculations, care must be taken to 
use realistic inputs for the wave rotor pressure ratio. 
Reference 1, p. 18, suggests some practical limits. As the 
static pressure ratio (SPR) increases, the efficiency of the 
wave rotor decreases. For good efficiency, total pressure 
ratios should be about 2.0. In the gas generator model, or 
pressure gain combustor, the total pressure ratio (TPR) should 
range from 1.1 to 2.0, corresponding to a total temperature 
(TTR) range between 1.5 to 3.8 [Refs 14; 15:p. 707]. 

The treatment of the wave rotor as an ideal component in 
the engine cycle is a limitation which eventually must be 
removed. The engine SFC and specific thrust depend greatly on 


the component efficiencies used in the performance 
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calculations. Their affects are not always obvious. A key 
parameter in turbofan engine design is the bypass ratio. As 
the bypass ratio increases, specific thrust decreases but SFC 
decreases to a minimum and then increases [Ref. 8:p. 163]. 
The optimum arrangement becomes dependent on the efficiencies 


attained by the individual components. 


B. RESULTS 

A study was conducted using WRCOMP and ENGINE to compare 
a turbofan with mixed exhausts and incorporating a wave rotor, 
with the baseline engine operating with the same turbine inlet 
temperatures. Table III lists the input parameters used in 
the comparison. A list of the component efficiencies is given 
in Appendix A. 

In maintaining a constant TIT, the effect of adding a wave 
rotor component to the turbofan was to change both the SFC and 
the specific thrust. The results obtained for an engine 
designed with a wave rotor component operating in a gas 
generator mode are shown in Figure 12. 

For each run the fan pressure ratio (n,;) was fixed and the 
compressor pressure ratio (n,.) was varied. Initially, the 
bypass ratio was allowed to vary to satisfy the equal-static- 
pressure requirement at the splitter plate, for a mixed 
exhaust. The calculated bypass ratio for all the test cases 
ranged broadly from zero to three. The overall pressure ratio 


was varied by fixing the fan and the wave rotor pressure 
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TABLE III 


RANGES OF INPUT PARAMETERS IN ENGINE COMPARISON 


0.79 BASELINE TURBOFAN WITH 
35000 ft TURBOFAN WAVE ROTOR 
1860 


Mier = 0. 
§ (bypass ratio) 


2.3 = 3.8 
ee ee 3.4 10.87 


Static Press. 
Ratio 





ratios and varying the compressor pressure ratio. The fan 
pressure ratio was then increased. At each fan pressure 
ratio, a specific range of values was obtained for the 
compressor pressure ratio for which valid solutions existed. 
As the fan pressure ratio was increased, the feasible range of 
values for the compressor pressure ratio decreased (the best 
performance range for the overall pressure ratio was from 22 
to 26), and the bypass ratio decreased. For example, for 
approximately the same TIT = 1860 and top, = 23.6 and for tn; = 
2.3 then nx, = 6.4 and § = 1.802; for nm, = 2.8 then n, = 5.3 and 
Bp - 1.144. With the wave rotor incorporated, the SFC was 


consistently lower than the baseline turbofan. The specific 
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Figure 12. Effect of Varying PRHC and PRLC for a 
Turbofan with a Wave Rotor 


(M, = 0.79, h = 35,000 ft, AMACH8 = 0.4, fB varies) 


thrust was also higher than the baseline engine except for the 
case Of nm, = 2.3. 
Since a large bypass ratio may not be acceptable for a 


cruise missile (a larger bypass ratio would require a larger 


fan size which may not meet the packaging requirements), the 





bypass ratio was set to B = 0.5, and the runs repeated. The 
results are given in Figure 13. A much smaller range of 


feasible compressor ratios resulted. 
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Figure 13. Effect of Varying Compressor Pressure 
Ratio at B = 0.5 for a Turbofan 


(M, = 0.79, h = 35,000 ft, AMACH8 = 0.4, f varies) 
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Table IV shows the trade-off between % change in SFC and 
specific thrust (referred to the baseline turbofan engine) as 


the fan pressure ratio was changed. 


TABLE IV 


PERCENT DIFFERENCE FROM BASELINE ENGINE 


+19.38 


. 
| se | = ao T9765 





Since it was not obvious that the improved performance was 
due to the wave rotor, the baseline engine performance was 
calculated as the fan pressure ratio was varied from 2.0 to 
3.5. The results are shown in Figure 14. It can be seen in 
Figure 14 that, as the fan pressure ratio was increased, the 
SFC increased very slightly and the specific thrust was also 
almost constant. Thus the improved performance in Table IV 


was a direct result of the wave rotor. 
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Figure 14. Effect of Varying the Fan Pressure Ratio 
for the Baseline Turbofan 


Table © summarizes the results of the comparison of the 
baseline turbofan with an engine that incorporates a wave 
rotor component. The baseline performance is shown compared 
to the wave rotor turbofan with the same turbine inlet 
temperature (TIT), the same bypass ratio (f), or the best 
specific thrust (ST). Of interest to a missile application, 
a 9% increase in specific thrust and a 7.3% decrease in SFC 


was obtained by incorporating a wave rotor into the baseline 
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turbofan and increasing the fan pressure ratio to 2.8, while 
maintaining $f = 0.5. 

Table V also infers other options available in particular 
applications. For example, if the capture area was not a 
design limitation, a higher bypass ratio (ff = 0.91) and fan 
pressure ratio (mn; = 3.0) could improve the specific thrust 
(ST) by 16.5% and decrease the SFC by 6.2% while operating at 
the same TIT = 1860. Similar improvements in specific thrust 
(ST) and SFC were obtained for lower values of TIT at the 


specified n;, as shown in the "best ST" entry in Table V. 
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TABLE V 


COMPARISON OF BASELINE TURBOFAN WITH WAVE ROTOR TURBOFAN 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


In the present study, the computer codes ENGINE and WRCOMP 
were installed and successfully operated on a VAX computer 
system. Modifications were made to enhance the utility of the 
programs and the results reported by the codes' author, A. 
Mathur, were reproduced. 

Comparisons between the predictions of ENGINE and ONX (by 
3. Mattingly) led to the conclusion that the use of constant 
specific heats in ONX leads to less accurate results than can 
be obtained with ENGINE, which includes real gas behavior 
throughout. Indeed, the predictions of ONX were found to be 
highly sensitive to the input values of specific heats. 

In first attempts to use the codes to examine the 
potential benefits of incorporating a wave rotor component in 
a mixed-exhaust by-pass fan engine, it was found that benefits 
were available which could not be obtained by varying the 
conventional engine components. For example, while 
maintaining the same by-pass ratio, a 9% increase in specific 
thrust and 7.3% decrease in SFC were predicted if, in addition 
to incorporating the wave rotor, the fan pressure ratio was 
increased from 2.3 to 2.8. Increases in specific thrust of 
almost 47% could be obtained with no change in SFC by 


increasing the fan pressure ratio to 3.8. 
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These results show that engines incorporating wave rotors 
can have improved performance over conventional gas turbines 
which are Jimited to the same turbine inlet temrersature 
levels. Further extensions of the ENGINE code 

(1) to include other engine types 


(2) to include friction and heat transfer in the wave rotor 
simulation 


are recommended. Also, studies to examine the effects of 
varying static pressure ratio and stagnation pressure ratio 
across the wave rotor, and of increasing turbine inlet 
temperature, should be carried out. 

Finally, an experimental program to validate the wave 
rotor flow predictions and performance levels used in cycle 


studies, is requirec. 
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APPENDIX A 


PROGRAM ENGINE 


PROCEDURES 
1. Log on to the VAX 
2. At the DCL prompt type: 
Sgksetup...ccsseces -.-.-initialize GRAFkit software 
Suis devon elevateneteS eioters ee sends plots to screen 
Srun ENGINEIA....... »...runs ENGINE program 
ore ehiotne. 3s vevensesioxe select options from screen 
3. To make changes in ENGINE1A.FOR: 
Sgksetup _.... ee oeee.---Only required once 
SULTS!) = swe bee ees a sears 
Sedit ENGINE1A. FOR ....to enter edit mode 
*C ip Siete aka .....to enter full page editing 
<CNTRL/Z> ...0sse0e2---0Xits edit mode 
*xexit oe eeeeeceeee -Saves Changes, DCL prompt returns 
SFORTRAN ENGINE1A .....compiles FORTRAN code 
Slink ENGINE1A,'GKL' ..links ENGINE1A and GRAFkit 
$run ENGINE1A 
see eeeceeeeeee-fOllow screen instructions 
4. To send your results to a laser printer: 


SGkSeGtup: stuccin initialize GRAFkit 
Wisi are 5oce as environment 
SINGS chee hea oes sends output to laser printer 
$define gk_out plot.dat .sends output to <fn>.<ft> 
Xe site ada wl acdae ale where plot.dat is your output 
ooeceeeeeeee efile name 
S1nN03S_ SP .oceesesece ....output to laser for square 
ce cece ceee plot 
$run ENGINELA seeeeeeeees-fOllow screen instructions 


Sprint/passall/que= 1n03 plot.dat 
See wee eee we -plot.dat is the same file above 


Note: Since the results are sent to an output file, 
"plot.dat", they will not appear on the screen. 
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B. POLYTROPIC EFFICIENCY SAMPLE CALCULATIONS 


The following sample calculations use the ENGINE output 
sample results from Section C in Appendix A. 
Compressor efficiency: > 


In terms of cycle parameters, 


ideal work interaction 


Ne ~ Sétual work interaction 


A(1) 
for a given pressure ratio. In terms of cycle parameters, 
Ved 


Mise oa A(2) 
tT -1 


Conpressor polytropic efficiency: 


For a given differential pressure ratio, 


ideal work interaction 


= 2S es ee A(3 
€ actual work interaction (3) 
The relationship with compressor efficiency is., 
You1 
Tome ote Te A(4) 
so that 
-1 
ins; Ye ar) ol i A(5) 
Ve : Pe 
or 
Y--1 int 
c= c Cc A(6) 
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For the fan (LPC): 


From ENGINE1.DAT: 


. = 1.4004 
q,= =o = 2.2941 
t,= 282 = 1.3108 

444 


giving, with Eq. A(6), 
e, = 0.8772 
For the high pressure compressor (HPC): 


From ENGINE1.DAT: 


Yo = 1.3870 
wee ree. = 10.8547 
ie eae 
giving, with Eq. A(6), 
e. = 0.9132 
Turbine efficiency: 
For a given pressure ratio, 
nee actual work interaction A(7) 
f ideal work interaction 


In terms of cycle parameters, 
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Turbine polytropic efficiency: 


For a given differential pressure ratio, 


actual work interaction 
ideal work interaction 


Qe. = 


The relationship to turbine efficiency is 





e, + (y¥_-1) 
2 Y 
tT, TT, : 
so that 
-1 
Int, =e, Yeo" an r, 
Yer 
or 
Int 
e.= Te oy t 


Y,71 in ®, 


For the high pressure turbine (HPT): 
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A(8) 


A(9) 


A(10) 


A(11) 


A(12) 





From ENGINE1.DAT: 





Y, = 1.3436 
22.8 

= = 0.1851 
t 123.2 
1295 

= = 0.6962 
f 1860 


giving, with Eq. A(12), 
e, = 0.8392 


For the low pressure turbine (LPT): 





¥_ = 1.3612 
11.6 

= ~ 0.5088 
t” 22.8 

_= 210 - 0.8525 
1295 


giving, with Eq. A(12), 
e, = 0.8899 
the polytropic efficiency for 


ENGINE was modified to calculate 


comparison with ONX. 
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SAMPLE RESULTS 


ae 


ON’ 12 xs 


ENCINE O14 


NOTATION COMPARISON FOR ONX AND ENCTWE 





Layout of baseline 


st 


Btion Identifications 


Turbofan 


O-1 Inlet Section 
1-2 Low Pressure Compressor 
2-3 High Pressure Compressor 
3-4 Mein Burner 
4-5 High Pressure Turbine 
8-6 Constant Pressure Mixer 
6-8 Low Pressure Turbine 
2-7 Bypass Duct 
7,8-9 Constant Area Mixer 
9-10 Afterburner 
10-12 Norzle 
0 3 4 5 6 7389 0 2 
Pf 


cer) 





Wave Rotor as Gas Generator 


Station IdentiFfiosatrions 


ch ve ve 


rue 
MQOONGAAA 


™ 
Lanvauaanna 


ne 

oe 
1 

ome 


Intet section 

Lov Pressure Compressor 
High Fressure Compressor 
Wave Potnr ond Main Burner 
High fresgure tendine 
Constant Pressure Hixer 
Low Pressure Turbine 
Bypess Duct 

Constant Area Hixer 
Afterburner 

Nozzle 
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BASELINE ENGINE CALCULATIONS 
a. DATA INPUT 
Flight parameters: 
M, = 0.79 
h = 35000 ft. 
P,.= 3-468 psi....\ (calculated within 


| the program) 


URE o-Ps i” eR 
Design Choices: 
Tre = 520 °R 
FLHV = 18000 BIU/1bm-k 


LPT exit Mach number (AMACH8) = 0 4 





Re 


Te 


PRLC = 2. 


PRHC 


10.8696 


Component Figures of Merit: 


Ne 
Ne 


Tepr = 
= ETALT 


N ipr 
Np 


ETALC 
ETAHC 
ETAHT 


ETACC 


tt oW 


ooo0o°o 


tt 


3 


. 86 
.88 
. 87 
-90 
96 


Maximum Component Temperatures: 


= 1860 


= 3400 


= 2085 
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°R 


°R 


°R 








Aircraft System Parameters: 


XI = HPT coolant air = 0.0 


Bypass ratio = B = 0.5 


Component losses: 


4 pp = BPLOSS = 0.02 
Ace = CCLOSS = 0.03 
ap = ABLOSS = 0.02 if AB off 


= 0.06 if AB on 


Miscellaneous: 
Timax = PIMAX = 0.97 
pmax = PIDMAX = 0.97 
VCOEFF = 0.97 

With the Wave Rotor: 

Tspp = SPR = 0.4 
Trp, = TPR = 1.6 
RREF= reference density in WR = 3.6 kg/m? 


Wr exit Mach number = AMEX = 0.9 


3.8 


It 
tt 


Te PRLC 


Te, = PRHC = 4.1 
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b. ONX OUTPUT 


HARES HHHKERREEE THERA EH 


MACH NO = 7790 

ALT (FT) = 35000. 

T CR) 2 394.10 

FO (FSIA) = 3.4680 

DENSITY 3.00073824 

{SLUG/CUFT) 

CF C = .243 BIU/LRM-R 

Cr YT = .2@65 BIU/LEM-K 

GaMNA C = 1.400 

GAMMA T = 1,350 

TT4 WAX = £860. CF) 

H - FUEL CEITU/LEM) = 18000, 

cio = .Qocdy 

COOLING AIR #1 = oO 4% 

COOLING AIR #2 = 200 4% 

FOsES = 1.00 

eee MIXER eae 
HEHEHE HER ENE EHRHREKERHHEER 

TAU R= 1.125 

FIR = 1,509 

FID = F770 

TAU L = S.138 

FUR TO oO KW 

TAU C’= 1.312 

EIA C= .8b04 

FIC = 95.90 

FIC’ = @,300 

TAU C*' = 1.312 

ETA C' = .B604 

FIOCH =10,870 

TAU CH = 2.198 

ETA CH = .6824% 

FIloTH = .1759 

TAU TH = . 6855 

ETA TH = .8649 

Fro’ = .63290 

TALE TO = 69014 

WITHO'UT AB 

F = .O116 

FO = .0116 

F/M = 42.415 LEF/LEM/S 

Ss = 49834 1/HR 

19/710 = 2.0535 

V9/VQ, = 2.7434 

MI/MDO = 1.9615 

AI/AN = 7480 

APF/AB = 1.28571 

THRUST = 6483. LEF 


ON--DESIGN CALCULATIONS 


TURBOFAN ENGINE WITH MIXED EXHAUST 


INFUT DATA 


RESULTS 
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AHKKHHHHERREPRHEREERHEEE 


AL FHA = 200 
FFIEC’ = 2.30 
FI D (HAX) = 297 
FIOR = 97 
FION = -97 
EFFICIENCY 

BURNCR = 95 
MECH HI FR = 1.00 
NFCH EQ FR = 1.00 
LF COMES (FAN) = 2-88 (EC) 
HF COME F. = »71 (CECH) 
HF TURRINE = »-84 (ETH) 
LF TURBINE = 89 CETL) 
FUR MECH EFF = 1,09 
FLEED AIR = 290 % 
FI MIXER MAX = 1.00 

RHEE REE EPAPER EMRE HEE 
AD = 973.1 FI/SEC 


Vei= 7468.7 FI/SEL 


MASG FLOW = 


eog,00 LEM /SEL 


AREA CERO = 19.954 SOFT 
APTA CERO# = 10,507 SOFT 
PTS /6O = 3,36? 

TIS /1TO = 1.476 

TAU ML = 1.0000 

TAU Mea = 1 ..000n 

TAU M = 1.0010 

FIOM = 1.0000 

Lia ba = 2 CIO 

NS 4 oO 

Ma = OOO 

AN SAS = ony 

GAMMA M= 1.20500 

CFM = 22659 

ETA TL = —B7923 

FIS/FO = 3.B494 


ENGINE OUTPUT 


FPR=2.3 OPR=25.0 TIT=1860. ABTEMP=3400. MACH NO=0.79 ALT=35000. 


STATION IDENT: ENGINE INLET, STATIONS 0-1 

Tle 394. TO- 444. PlIe 3.5 POw 5.1 FI|N0.000 FO=.000 vo= 769. 
STATION IDENT: LP COMPRESSOR, STATIONS 1-2 

Tile 444. TO= 582. PIe 5.1 PO= 11.7 Fre.000 FO-.000 
ETALC#0.8600 GBAR=1.4004 

ELPC*0.6757 

STAION IDENT: HP COMPRESSOR, STATIONS 2-3 

Tile 582. TO-1206. Pile 11.7 PO#127.0 Fie.000 FO-.000 
ETAHC=0.8800 GBAR=1.3670 

EHPC=0.9144 

STATION IDENT: MAIN BURNER, STATIONS 3-4 

T121206. TO-1860. PI|el127.0 PO9123.2 FIe.000 FOe.010 
STATION IDENT: HP TURBINE, STATIONS 4-5 

T1e1860. TO-1295. Pre|el123.2 PO= 22.8 FI=.010 FO=-.010 
ETA=0.8700 GDBAR@1.3436 

ET=0.8379 

STATION IDENT: PRE LPT MIXING, STATIONS 3,5-6 

The1206. 7221295. TO=1295. PO» 22.8 Fil=.000 F2=.010 FO=.010 
STATION IDENT: BYPASS DUCT, STATIONS 2-7 

Tle 582. TO= 582. PYI= 11.7 PO= 11.5 FI=.000 FO=.000 
STATION IDENY: LP TURBINE, STATIONS 6-8 

TI#1295. TC#'104. PIs 22.8 PO= 11.6 Fie.010 FO=.010 
ETA-0.9000 GEAR=1.3612 

eT-0.8904 

STATION IDENT: CONST. AREA MIXER, STATIONS 7,8-9 

T1e1104. T2= 582. TO 939. Ple 11.6 P2= 11.5 PO= 11.2 
Fie.010 F2=.000 FO=.007 M1"0.40 M2=—0.37 MO=-0.40 ARAT=0.38 
STATION IDENT: AFTERBURNER, STATIONS 9-10 

Tre 939 TO 939. Pre 11.2 Pow 11.0 FIe.007 FO@#.007 
STATION IDENT: EXHAUST NOZZLE, STATIONS 11-12 


Tle 939. TO= 697, PIe= 11.0 POe 3.5 FIe.007 VO01l728. MOW1.34 
GLOCO#1.3948 PRN]=.31661 AEXIT#1294.1 


PERFORMANCE PARAMETERS 


AMACH8e 0.400 


SP.THRUST = 39.18 SFC = 0.832 BETA= 0.500 EXMACHe 1.336 
(LBF/LBM/S) (LBM/HR/LBF ) 
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3. TURBOFAN WITH WAVE ROTOR: ENGINE OUTPUT 


FPR=3.8 OPR=-24.9 TIT#=1860. ABTEMP=3400. MACHR NO=0.79 ALT=35000. 
STATION IDENT: ENGINE INLET, STATIONS 0-1 
Tie 394. TOe 444. Pre 3.5 pom §.1 FI-0.000 FO=§.000 vom 769. 
STATION IDENT: LP COMPRESSOR, STATIONS 1-2 
Tle 444. TO= 682. Pre 5.1 poe 19.3 FI=.000 Fo-.000 
ETALC=0.8600 GBAR=1.3992 
ELPC=0.88641 
STAION IDENT: HP COMPRESSOR, STATIONS 2-3 
Tie 682. TO=1056. Pre 19.3 Pow 79.2 FI=.000 FO=8.000 
ETAHC=0.6600 GBAR=1.3870 
EWPC=0.9022 
STATION IDENT: WAVE ROTOR, STATIONS 3-4 
TCCIN~ 1199.0 TO= 1972.7  TREF= 2197.7 
PIN= 79.2 PO= 122.9 PREFe 183.1 
FIN=0.0000 FO=0.0165 GLOCOe1.371 IDBi#1 
STATION IDENT: HP TURBINE, STATIONS 6-5 
TI"1973. TO:1649. P1#122.9 pow 52.6 Fi=.016 FO=.016 
ETA=0.8700 GBAR#=1.3263 
£T=0.8551 
STATION IDENT: PRE LPT MIXING, STATIONS 3,5-6 
T1e1056. T2~1649. TO—1649. PO= 52.4 Fil=.000 F2~.016 FO=.016 
STATION IDENT: BYPASS DUCT, STATIONS 2-7 
Tle 662. TO= 682. PI= 19.3 PO= 18.9 Fre.000 FO-.000 
STATION IDENT: LP TURBINE, STATIONS 6-8 
TI©1649. TO|91335, Pile 52.4 PO= 20.6 Fie.016 FO#.016 
ETA#0.9000 GBAR@1.3416 
ET=0.8866 
STATION IDENT: CONST. AR MIXE%, STATIONS 71,86-9 
T1#1335. T2= 662. TO-1133. Pie» 20.6 P2= 18.9 Po= 19.4 
Fl=.016 F2=.000 FO~.011 m1=0.40 M2e0.19 MO~0.31 ARAT=0.74 
STATION IDENT: AFTERBURNER, STATIONS 9-10 
T1e1133. 7TO91133. Pre 19.4 POH 19.0 Fie.011 FO=.011 
STATION IDENT: EXHAUST NOZZLE, STATIONS 11-12 


Tre1133, TO= 737. PIe 19.0 Poe 3.5 FIe.011 VO—02232. mO=1.68 
GLOCO=1,3690 PRN=.18289 AEXIT#1329.8 


PERFORMANCE PARAMETERS 


AMACHB= 0.400 
SP.THRUST = 46.22 SFC = 0.856 BETA 0.500 EXMACHe 1.678 


{LBF/LBM/S) (LBM/HR/LBF ) 
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D. 


ENGINE 


C fan 
c 


11 


PROGRAM LISTING 


PROGRAM ENGINE 
DIMENSION XARAY(125),YARAY(125), ZARAY(125)},B8PR(125) 
REAL XP1(50),YP1(50),YP2(50) 
REAL ST,SPTHR,SFC,PRAT 
INTEGER CONFIG, MCONFIG,II,33,11,01 
CHARACTER*6 XiILABEL 
CHARACTER*15 YILABEL 
CHARACTER*3 Y2LABEL 
COMMON/ETAS/ETALC, ETAHC, ETACC, ETAHT, ETALT, ETAAB, VCOEFF 
COMMON/LOSSES/DIFLOS,CCLOSS,ABLOSS, BPLOSS 
COMMON/BURN/HFUEL, FLAY, TMAXCC, TMAXAS, THAXWR 
COMMON/WORK1 /WORKLC, WORKHC 
COMMON/CONST/CJ,G 
COMMON CONFIG 
DEFINITION OF PROGRAM SWITCHES 
CONFIG=O0:CONVENTIONAL ENGINE; CONFIG=1:WAVE-ROTOR ENGINE. 
KODE=0:NON-AFTERBURNING CASE; RODE=1: AFTERBURNER LIT CASE. 
KASET=0:B8YPA8S RATIO (BETA) PRESCRIBED. KASET@1:BETA CALCULATED. 
DATA 1DC,1DB,IDT,IDM/0,0,0,0/ 
DATA XP1/50*0.0/ 
DATA YP1,50*0.0/ 
DATA YP2/50*0.0/ 
IFLAG=-0 
INPUT DATA (GENERAL) 
RERRARAERRARAATEHARHARAARAHRERARRARAARRAERRARAHRARKRARARAARRARARAEAHOAE 
INITIALIZE THF DATA BASE s**t ata teat tA Rh kAHAAHAHHARHEAARAEAEKAEhEh 
WRITE(6,*)’WOULD YOU LIKE TO SEE THE POSSIBLE CONFIGURATIONS 
> FOR ENGINE!? (l=YES ONO)’ 
WRITE(6,*)'Hit <Return> for next picture and to continue. 
READ(S5,°(1I)* )IDESIGN 
IF (IDESIC™*.£Q.0) GOTO 11 
CALL LIBSSPAWN ('RENDER £NG1.UIS’) 
CALL LIBSSPAWN (/ RENDER ENG2.UI8" ) 
CONTINUE 
WRITE(6,*)°INPUT DESIRED CONFIGURATION:’ 
WR1TE(6,*)'CONFIG=(OeCONVENTIONAL L=WAVE ROTOR)’ 
READ(5,‘'(1)'’ )CONFIG 
MTO"FIG*CONFIG 
WRITE(6,*) "2ZODE@{0=NON A/B 1<A/B ON)’ 
READ(5,°(1)") ROLE 
WRITE(6,*)'RASET=(O=sBYPASS RATIO PRESCRIBED leBETA CALCULATED)’ 
READ(5,'(1)") KASET 
WRITE(6,*)*WOULD YOU LIKE TO USE YOUR OWN PARAMETERS OR THE 
> DEFAULT’ 
WRITE(6,*)’ PARAMETERS? (1=OWN, O@DEFAULT)’ 
READ(5,'(1)’)1Q1 
Ir (1Q01.£Q.1) THEN 
CALL FLIGHT(FMACH,ALT,PSO,TSO) 
WRITE(6,*)°CONFIG=" , CONFIG 
If (70"N"'G. 869.1) THEN 
CALL WRSPECS(S5PR, TPR, RREF , AMEX) 
ENDIF 
WRITE(6,*)*CONFIG=",CONFIG 
CALL DESLIM(TFUEL, FLAV,AMACHS, HFUEL, PRLC, PRAC, PRAT, TPR, 
>CONFIG) 
CALL COMPFOM(ETALC, ETAHC, ETAHT, ETALT, ETACC, ETAAB) 
CALL TMAX(TMAXCC, TMAXAB, TMAXWR) 
CALL ACSYS(CTO,TMAXCC,X1,BETA) 
CALL COMLCSS(BPLOSs,CCLOSS,ABLOSS) 
CALL COMPETA( PIMAX, PIDMAX, VCOEFF) 
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ELSEIF (1Q1.£9.0) THEN 
CALL BASE1(FMACH,ALT,PSO,TSO, TFUEL, FLAV,AMACHB,AFUEL, 
>PRLC, PRHC, PRAT, ETALC, ETAHC, ETAHT, ETALT, £TACC, ETAAB, EATDB, 
>TMAXCC, THAXAB, TMAXWR, CTO,X1,BETA, BPLOSS , CCLOSS, ABLOSS, 
>PIMAX, PIDMAX, VCOEFF, SPR, TPR, RREF ,AMEX, CONFIG, KODE} 
ENDIF 
c~---- SET INTIAL CONSTANTS8-----~--~----~--------~-------~~ a oe 


Cee eee 668888888 8S FOE SENSO SOS ESOS SERS SSA AH SEAS ASSSeeesnaseasess 
C----- END OF INPUT DATA 
Cees ee ek 6886 68 ESOS SOE SOS SENSES SESE SORES BESS SS SSS SH SASS SSeer Tee Ssee8e 
OPEN (UNIT#=9,FILE@="’ENGINE1.DAT’ ,STATUS="NEW’ ) 
Cowen weseaeenene saan eas sasee ae Seen nee un een asaceresecusseanneneenauaae 
C *** CALLING SEQUENCE OF COMPONENT SUBROUTINES FOR CHOSEN ENGINE *#* 
Cornre CONFIGURATION 
Cee 88 OS 6 OS CEN O ENS A USES S BSSSCHSSEETAHRSSS VOUS SAT Ses SASS UENES ABSCESS 
WRITE(6,*)"WOULD YOU LIKE TO ITERATE A VARIABLE (1=YE6 O=NO)’ 
READ(5,’(1)')IANS1 
IP (IANS1.£Q.0) THEN 


GO TO 6 
ELSE 

WRITESE SOPOT THE DESIRED ITEPATION PARAMETER (ie 1)’ 
WRITE(6,%)’ 1. BYPASS RATIO’ 

WRITE(6,*)’ 2. FAN PRESSURE RATIO’ 

WRITE(6,*)° 3. COMPRESSOR RATIO’ 

WRITE(6,%)’ 4. LPT EXIT MACH NUMBER’ 

READ(5,*)IANS4 

IANSW4eIANS4 


WRITE(6,*)'INPUT MIN VALUE =' 
READ(5,*)S8MIN 
WRITE(6,*)'" INPUT MAX VALUE =’ 
READ(5,*)SMAX 
WRITE(6,*)" INPUT STEP SIZE =’ 
READ(S,*)STEP 
OPEN(UNIT=14,FILE='’ PLOT1.DAT’ ,STATUS=' NEW’ ) 
TF (IANS4.£Q.1) THEN 
XILABEL©’ BETA ’ 
¢ WRITE(14,°(2X,A,2X,A,2X,A)')* BETA’,’ 8T ',° BFC’ 
ELSEIF (IANS4.£EQ.2) THEN 
XILABEL©=' PRLC ' 
c WRITE(14,°(2X,A,2X,A,2X,A)")’ PRLC’,’ ST °,% SFC ’ 
ELSEIF (IANS4.£Q.3) THEN 
XLILABEL@’ PRHC ° 


C WRITE(14,°(2X,A,2X,A,2X,A)°)’ PRHC’,’ 8ST °," SFC ' 
ELSEIF (IANS4.£0.4) THEN 
XILABEL@='AMACHB’ 
¢c WRITE(14,°€2X,A,2X,A,2X,A)")'"AMACHS’," 8T °," BFC ' 
ENOIF 


WRITE(14,°(2%,R,2X,A,2K,A,2K,A,2K,A,2X,A6,2%,A,2K,A)')’ BETA ’, 
>’ PRLC ',’ PRHC °,’ PRAT ',’° TIT ',XILABEL,*° 8T °,’ BFC ’ 
1C2e1 
C------START LARGE LOOP ITERATION FOR ITERATION VARIABLE~-----~--- 
DO 41 SVALUE“SMIN, SMAX,E6TEP 
IF (IANS4,£Q.1) THEN 
BETA*SVALUE 
ELSEIF (IANS4.E0.2) THEN 
PRLC@SVALUE 
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ELSEIF (IANS4.£Q.3) THEN 


PRHC@SVALUE 
ELSEIF (1ANS4,£Q.4) THEN 
AMACH6=SVALUE - 
ENDIF 
IF (CONFIG.£Q.1) THEN 
PRAT=PRLC* PRHC*TPR 
ELSE 
PRAT@PRLC* PRAC 
ENDIF 
C----- WRITE STATION RESULTS TO DEVICE 9 ON FILE ENGINE). DAT------------ 
C----- ITERATION PLOT RESULTS ARE SENT TO 14 AND PLOT1.DAT-------------- 


6 WRITE(9,1)PRLC, PRAT, TMAXCC, TMAXAB, FMACH,ALT 
1 FORMAT(2X,°FPR='’,F3.1,2X,'OPRo',F4.1,2%, 'TIT=’,F5.0,2K%, "ABTEMP~’, 
>F5.0,2X,'’MACH NO=’,F4.2,2X,'ALT=',F6.0/) 
CALL INLET(PSO,TSO,FO,P1,T1,F1,VO, FMACH, PIDMAX) 
CALL COMPR(PRLC,P1,T1,P2,T2,T218,F1,F2,IDC) 
CALL COMPR(PRHC,P2,T2,P3,T3,T3I8,F2,F3, IDC) 
IF(CONFIG.£Q.1) GoTo 10 
CALL BURNER(P3,T3,P4,T4,F3,"F4,1DB, KODE, CONFIG) | 





._ GOTO 20 
ce 10 CALL ROTOR(P3,P4,T3,T4,PREF,RREF,TREF,F3,F4,SPR,TPR, AMEX, TCCIN, 
¢ >IDB, CONFIG) 


10 CALL ROTOR(PRAT,71,P3,P4,T3,74,PREF,RREF, TPEP, F3,F4,5PR,TPR, AMEX, 
>TCCIN, IDB, CONFIG) 
20 CALL TURBIN(P4,T4,P5,T5,T518,F4,F5,85,PRAT,PRLC, PRAT, PRLT,XI,BETA, 
>RASET,IDT) 
CALL MIXCP(F3,F5,P3,73,P5,T5,86,P6,T5,K1,IDM) 
CALL DUCT(P2,T2,F2,P7,T7,F7,H7) 
CALL TURBIN(P6,T6,P8,T8,T81S,F6,F8,H8,PRAT,PRLC, PRAT, PRLT,XI, BETA, 
>RASET,1IDT) 
IF(KASET.£Q.1) GOTO 30 
CALL MIXCA(P7,P8,T7,T8,H7,H8,f7,F8,BETA,AMACHS, PIMAX,T9,£9,P9, 
>AMACHY) 
GOTO 40 . 
30 CALL MIXCP(f7,F8,P7,T7,P8,T8,F9,P9,T9,BETA, IDM) 
40 CALL BURNER(PS,T9,P10,T10,F9,F30,1D8, KODE, CONFIG) 
CALL NOZZLE(P10,T10,P12,712,T1218,F10,PRLC, VEXIT, EXMACH, P50, IFLAG) 
IF(IFLAG.EQ.1) GOTO 50 
CALL PERF (AMACH8,F10,VO,VEXIT, BETA, SPTHR, SFC, EXMACH) - 
WRITE(14,3014) BETA, PRLC,PRHC, PRAT, T4,SVALUE,SPTHR,SFC 
WRITE(14,°(2X,F6.3,2X,F6.3,2X,F6.3)" )SVALUE, SPTHR, SFC 
WRITE(6,*)'IC2=',1C2 
WRITE(6,*)’ST=’,SPTHR, SFC=’,SFC 
XP1(1C2) = SVALUE 
YP1(IC2) = SPTHR 
YP2(1C2) = 8FC 
c WRITE(6,*)'XPl=’,XP1(IC2),'¥Pl=’,YPI(IC2),'¥P2@’ ,¥P2(1C2) 
If (SVALUE.LE.SMAX) 1C2=1C241 
t--~--------- RESET THE COUNTER SWITCHES EACH LOOP 
rpc=0 
IDB=0 
IDT=0 * 
IDM=0 
CONFIG=MCONFIG 


aad 


c WRITE(6,*)°IDC=',IDC,’IDB=',IDB,’IDT=’, IDT, "IDM=’, IDM 
c WRITE(6,*)’ CONFIGe’ ,CONFIG 
41 CONTINUE 
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END IF 
ENDFILE 14 
REWIND 14 
CLOSE (UNIT=14) 
C8 8888 6 Ok 66 8 6 OES SEES E8888 SSO SSS SESS SOS ESSE Sse OE SHS SeSeSS 
c CALL GRAFKIT PLOT ROUTINE 
C88 88 888868 8888 8S CEES E SOO O8 STO ESS SOS SOS SSeS ES Sees eee aeEesese 
WRITE(6,*)’WOULD YOU LIKE TO PLOT THE RESULTS? (1=YES O=NO)’ 
READ(5S,°(1)')1C5 
IF (1C5.£9.1) TREN 


C----- NOW CALL GRAPHICS ROUTINE 
CALL GRSPLOT1(XP1,YP1,YP2, IANSW4,1C2-1)} 
ENDIF 
Cc aan ARAAARAARARARRRABREARAARARARARRAREARARRARRARAKEHARARHARAARHERARRERAR 
C *#* ***END OF USER-DEFINED CALLING SEQUENCE*### Ss tenannannnnaanane 
C ane RARARRHRARAERARRARARAARAAARARAREEARKRAREARERHARARARARARKARARRERER 


OPEN(UNIT=10,FILE*' ENGINE2.DAT’ ,STATUS='NEW’ ) 
WRITE(10,3000)P6,P7,P8,Pi1 
CLOSE(UNIT#=10) 
3000 gt ONAL ER ROS REDS AB Ble a Ee he Peta sete s Aire tie’ 4 
>F5.1//) 
3014 FORMAT(2X,76.3,2X,F6.3,2%,F6.3,2X,F76.3,2X,F6.1,2X%,F6.3,2X, 
>F6.3,2X%,F6.3) 
GOTO 70 
50 WRITE(6,60) 
60 FORMAT(5X, NOZZLE PRESSURE RATIO I8 NEGATIVE’ /) 
70 CLOSE(UNIT#=9) 
STOP 
END 
C88 06 6 O66 8 660988 SOS OOS SESE SEE ENOTES SS SEAS OTST SSeS SSS SeseenseVesenESs 
SUBROUTINE BASE1(FMACH,ALT,PSO,TSO, TFUEL, FLHV,AMACHB,HFUEL, 
>PRLC,PRHC, PRAT, ETALC, ETAHC, ETAHT, ETALT, ETACC, ETAAB, EATDB, 
>TMAXCC, Thanr.”, sMAYYE, CTO,XI,BETA, oP LOSS,CCLOSS,ABLOSS, 
>PIMAX, PIDMAX, VCOEFF, SPR, TPR, RREF ,AMEX, ICONEIG, IKODE) 
COMMON CONFIG 
C----- FLIGHT PARAMETERS 
FMACH#0.79 
ALT=35000. 
PSO=3.468 
TSO“394.1 
C----- WAVE ROTOR SPECS 
SPR-0.40 
TPR=1.60 
RREF=3.6 
AMEX=0.9 
C--~------ TCCIN IS SOLVED FOR IN ROTOR ROUTINE---------------- 
c TCCIN@1213.3 
Cw---- DESIGN/FUEL CHOICES 
TFUEL=520. 
FLHV=18000. 
AMACHB@0.4 
BFUEL@TFUEL*0.5 
If (ICONFIG.EQ.1} TREN -. 
PRLC#=3.8 
PRHCe4,1 
ELSE 
PRLC#2,3 
PRHC#10.8696 
ENDIF 
IF (ICONPIG.£Q.1) TREN 
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PRAToPRLC*PRAC* TPR 
c WRITE(6,*)'PRAT@=’,PRAT,'TPR=’, TPR,’ CONFIG@=!, ICONFIG 
GLseE 
PRAT=PRLC* PRHC 
ENDIF 
Con--- COMPONENT FOM 
ETALC=0. 66 
ETAHC=0.88 
ETAHT=0.67 
ETALT<-0.90 
ETACC#0.96 
IF (IKODE.£Q.0) THEN 
ETAAB=1.00 
ELSE 
ETAAB=0.96 
ENDIF 
C----- MAX COMPONENT TEMPERATURES (R) 
¢c NOTE: TMAXCC AND TMAXWR GET RESET TO VALUES CORRESPONDING TO 
c THE WAVE ROTOR ANALYSIS : TMAXCC@TREF, TMAXWReTREF 
TMAXCC#1860,0 
TMAXAB=3400. 
TMAXWR=2065 
Coonnn. AIRCRAFT SYSTEM PARAMETERS 
cTo-0.0 
WRITE(6,*)’ INPUT & COOLING OR BASE IT ON TMAXCC:’ 
WRITE(6,°%)’ (1=sINPUT FOR XI O=-BASE XI ON TMAXCC)’ 
READ(5,*)ICOOLI 
IF (ICCOL1.2£Q.0) THEN 
iF (TMAXCC.GT.240U) THEN 
K1I=2.C*( TMAXCC-2400.0)716000.0 
ELSE 
X1=0.0 
ENDIF 
ELSE 
WRITE(6,**’ INPUT XI’ 
READ(5,*)XI 
ENDIF 
BETA#0.5 
Co---- COMPONENT LOSSES 
BPLOSS=0.02 
CCLOSS#0.03 
2? (IRODE.EQ.0) THEN 
ABLOSS=0.02 
ELSE 
ABLOSS=0.06 
ENDIF 
Cxee== COMPONENT EFFICIENCIES 
PIMAX=0.97 
PIDMAX=0.97 
VCOE@F=0,4. 
RETURN 
end 
Ce 86868 006 6888 08 89S 680988285 OSCE SEER SO COTS RSS SEES SU ES SAE He CeEeseeRse 


c INPUT DATA BASE SUBROUTINES easecceeneuecn 


SUBROUTINE FLIGRT(FMACH,ALT, PSO, TSO) 

DIMENSION TT(0150),PP(0:50) 

DATA TT/1.0000,0.9931,0.9863,0.9794,0.9725,0.9656,0.9583,0.9519, 
20.9540,0.9381,0.9313,0.9244,0.9175,0.9107,0.9038,0.8969,0.8901, 
>0.6632,0.8764,0.8695,0.8626,0.8558,0.8489,0.86420,0.8352,0.6233, 











>0.6215,0.8146,0.8077,0.8009,0.7940,0.7872,0.7803,0.7735,0.7666, 
>0.7598,0.7529,0.7519,0.7519,0.7519,0.7519,0.7519,0.7519,0.7519, 
i >0.7519,0.7519,0.7519,0.7519,0.7519,0.7519,0.7519/ 

DATA PP/1.0000,0.9644,0.9298,0.8963,0.8637,0.8321,0.8014,0.7717, 
>0.7429,0.7149,0.6878,0.6616,0.6362,0.6115,0.5877,0.5646,0.5422, 
>0.5206,0.4997,0.4795,0.4599,0.4410,0.4227,0.4051,0.3880,0.3716, 
>0.3557,0.3404,0.3256,0.3113,0.2975,0.2843,0.2715,0.2592,0.2474, 

zs >0.2360,0.2250,0.2145,0.20464,0.1949,0.1858,0.1771,0.1688,0.1609, 
>0.1534,0.1462,0.1046,0.1329,0.1267,0.1208,0.1151/ 
WRITE(6,*)°INPUT FLIGHT CONDITIONS’ 
WRITE(6,*)’ INPUT FLIGHT MACH NUMBER: ’ 
READ(5,*) FMACH 
WRITE(6,%)° INPUT ALTITUDE (0-50000 FT)’ 
READ(5,*)ALT 
IALTeALT/1000 
TS0-(459.67+59.0)*TT( IALT) 
PS0014.696*PP(IALT) 
RETURN 
END 


SUBROUTINE DESLIM(TFUEL, FLAV,AMACHGA, HFUEL,PRLC,PRHC, PRAT, TPR, 
>ICONFIG) 

COMMON CONFIG 

WRITE(6,*)’ INPUT DESIGN/FUEL LIMITATIONS!’ 

WRITE(6,*)° TFUELe(TEMP OF FUEL IN R)’ 
READ(S,*) TFUEL 

WRITE(6,%)° FLAV=(FUEL LOWER HEATING VALUE)’ 
READ(5,*)PLAV 

WRITE(6,*)° AMACH6=(EXIT MACH AT LPT)’ 
READ(5,*)AMACHSA 

WRITEIE #)F PRLCH’ 
READ(5,*)PRLC 

WRITE(6,*)! PRHCa’ 
READ(5,*)PRHC 
IF (ICONFIG.£EQ.0) THEN 

TPR#1.0 

ENDIF 

PRAT=PRHC*PRLC* TPR 

HFUEL@TFUEL*O.5 

RETURN 
END 


SUBROUTINE COMPFOM(ETALC,ZTAHC,ETAHT,ETALT,ETACC,ETAAB! 
- WRITE(6,*%)’INPU1 COMPONENT FIGURES OF MERIT:(1=OWN O#STD)’ 
WRITE(6,*%)'(STO: ETALC=0.66 ETAHC-0.88 ETAHT=0.67 ETALT=0.90 
> ETACC#0.96 ETAAB=0.96)’ 
READ(S,°(1)’ )ISUB2 
IF (18UB2,£0.0) THEN 


ETALC=0.66 

ETAHC#=0.68 

ETAHT=0.67 

ETALT=0.90 

ETACC#0.96 bs 
ETAAB~0.96 

ELSE 

WRITE(6,#)° INPUT ETALC=’ 
READ(S,* )ETALC 

WRITE(6,*)? INPUT ETARCe’ 
READ(5,*)ETAHC 

WRITE(6,#)? INPUT ETAHT@=" 
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READ(5,*)&TART 
WRITE(6,*)’ INPUT ETALT#=’ 
READ(5,*)ETALT 
WRITE(6,*)° INPUT ETACC@' 
READ(5,*®)EBTACC 
WRITE(6,*)' INPUT ETAAB}’ 
READ(5,*)BTAAB 
enolr 
RETURN 
END 
c 
BUBROUTINE TMAX(TMAXCC, TMAXAB, TMAXWR ) 
WRITE(6,*) ‘INPUT MAX COMPONENT TEMPERATURES’ 
WRITE(6,*)* | INPUT TMAXCC=' 
WRITE(6,*)° TMAXCC 16 RESET TO TMAXWR IN ROTOR CONFIG’ 
READ(5,*)TMAXCC 
WRITE(6,*)’ | INPUT TMAXAB=° 
READ(5, *) TMAXAB 
WRITE(6,*)’ INPUT TMAXWR=’ 
WRITE(6,4)° " TMAXWR IS RESET TO TREF IN ROTOR CONFIG’ 
READ(S, #) TMAXWR 
RETURN 
END 
¢c 
SUBROUTINE ACSYS(CTO,TMAXCC,XI,BETA) 
WRITE(6,*%)'VERICLE SYSTEM PARAMETERS’ 
c WRITE(6,*)? INPUT CTO= (TAKE-OFF POWER)’ 
c READ(5,*)CTO 
WRITE(6,*)' | INPUT %& COOLING OR BASE IT ON TMAXCC:’ 
WRITE(6,¢)* (1sINPUT FOR XI O~BASE XI ON TMAXCC)’ 
READ(5,*)ICOOL1 
IF (I1COOL1.£0.6) THEN 
IF (TMAXCC.GT.2400) THEN 
X1-2.0*( TMAXCC-2400.0)/16000.0 
ELSE 
X1=0.0 
ENDIF 
ELSE 
WRITE(6,*) "INPUT XIe@! 
READ(S,*)XI 
ENDIF 
WRITE(6,*)’ INPUT BETA= (BY-PASS RATIO)’ 
READ(5,*)BETA 
RETURN 
END 
c 


SUBROUTINE COMLOSS(BPLOSS,CCLOSS,ABLOSS) 
WRITE(&,*)’COMPONENT LOSSES (PERCENTS):(1=OWN O=STD)’ 
: WRITE(6,*)'STD: BPLOSS=0.02 CCLOSS=0.03 ABLOSS=0.02 )’ 
READ(5,°(1}")105 
IF (105.£Q.0) TREN 
BFLOSS=0.02 
CCLOSS=0.03 ° 
ABLOGB=0.02 
ELSE 
WRITE(6,*)* INPUT BPLOSS=' 
READ(S,*)BPLOSS 
WRITE(6,*)* | INPUT CCLOSS~’ 
READ(S,*}CCLOSS 
WRITE(6,*)' INPUT ABLOSBe’ 
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ec 





READ(5,*)ABLOSS 
ENDIF 


RETURN 
END 


SUBROUTINE COMPETA( PIMAX, PIDMAX, VCOEFF) 


WRITE(6,*)'MISCELLANEOUS COMPONENT CHOICES: (1=OWN OSTD)’* 
WRITE(6,*)°(STO: PIMAX=0.97 EIDMAX=0.97 VCOEFF#0.97)" 
READ(S,°(1)'°)1Q6 
IF (106.£0.0) THEN 

PIMAK=0.97 

PIDMAX=-0.97 

vCOEFF=-0.97 


INPUT PIMAX=' 
IMAX 

INPUT PIDMAX=’ 
IDMAX 

INPUT VCOEFF= (NOZZLE EFF)’ 
COEFF 


<s< se 3 


BUBROUTINE WRSPECS(SPR,TPR, RREF, AMEX) 


an 


’ 
Conon- -TO INPUT T 
« 


WRITE(6,*)’NAVE ROTOR SPECIFICATIONS’ 
WRITE(6,*)’ INPUT SPR’ 


’ INPUT TPRe’ 
INPUT RREF=  kg/a‘3' 


WRITE(6, 
READ(S, 


INPUT TCCINe’ 
TCCIN 


RETURN 
END 


c 
C ete 


SHRARAHRARARARERRERAAARARRRARAEARARERHRAAAEARHAHERARARHHEERERARHAE 


SUBROUTINE INLET(PIN, TIN, FIN, PO,TO,FO,VO, FMACH, PIDMAX) 
COMMON/LOSSES/DIFLOS,CCLOSS,ABLO5S,BPLOSS 
COMMON/CONST/CJ,G 

RINSENTHFN( FIN, TIN) 

PHIIN@ENTREN(FIN, TIN) 


GLOCIN=GLOC(FIN, TIN, HIN) 
PTIN@-PIN®*(1.4(GLOCIN-1.)/2.*FMACH**2,)**#(CLOCIN/(GLOCIN-1.)) 
TTIN@TIN®(2.4(GLOCIN-1,)/2.*FMACH**2. } 

PID@PIDMAX 

IF(FMACH.GE.1.0) PID@PIDMAX*(1.-0.075*( FMACH-1.)*#1.35) 
DIFLOS#1.-PID 

PO-PTIN®*(1-DIFLOS) : S 

TO-TTIN 

ro=0.0 


BO-HIN 

PHIO=PRIIN 
VOeFMACH*SQRT(GLOCIN* RGAS( FIN) *TIN*G) 
WRITE(9,10) 
WRITE(9,20)TIN,TO,PIN, PO, FIN, FO,VO 
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20 FORMAT(5X,’STATION IDENT: ENGINE INLET, STATIONS 0-1°/) 
20 FORMAT(5X,'TI=',F5.0,2X,'TO=',F5.0,2X,' Plo’, F5.1,2K,'PO=", 
>P5.1,2X,'FI=’,5.3,2X,"FO=’, F4.3,2%,'VOm!,F5.0/) 
RETURN : 
END 
¢ BAR ARHARAARARARAHRAERARERAARARAHRARERHRERARARARERHRARARARHRARAARARAREERRAREERE 
SUBROUTINE COMPR(PR,PIN, TIN, PO, TO, TOIS,fIN,¥O,I0C) 
COMMON/ETAS /ETALC, ETARC, ETACC, ETAHT, ETALT, ETAAB , VCOEFF 
COMMON/LOSSES/DIFLOS, CCLOS8, ABLOSS , BPLOSS 
COMMON/WORK1/WORKLC, WORKEC 
COMMON/CONST/CJ,G 
IDC=1pC+1 
IF(IDC.£Q.1) GO TO 30 
C AOeeeeeCOMPUTE APC 
HIN@ENTHFN( FIN, TIN) 
PHIIN@ENTREN( FIN, TIN) 
PRIOIS=PHIIN+RGAS( FIN) /CJ*ALOG(PR) 
CALL PHIFN(FIN, PHIOIS, TOIS) 
HOIS@ENTHFN( FIN, TOIS) 
RO-HIN+(HOXS-HIN) /ETAHC 
WORKHC#RO-RIN 
CALL AFN(FIN,#O,TO} 
PHIO*ENTREN( FIN, TO) 
TO<FIN 
PO=PIN*PR 
GBAR=GAVG( FO, TIN, TOIS, PRIIN, PRIOIS) 
TEMPHC1*ALOG(TO/TIN) 
TEMPHC2=ALOG( PO/PIN) 
EHPC=(GBAR-1.0)/GBAR* TEMPHC2/TEMPHC1 
WRITE(6,*)'TO=', TO, TIN]©’, TIN, 'PO=’, PO,’ PIN@’, 
>PIN, ETAHCe’, ETARC, 'GBAR=’ ,GBAR 
WRITE(6,*)'EHPC=}", BHPC 
WRITE(9,10) 
WRITE(9,20)TIN, TO, PIN, PO, FIN, FO, ETAHC, GBAR, EHPC 
10 FORMAT(5X,’STAION IDENT: RP COMPRESSOR, STATIONS 2-3'/) 
20 =FORMAT(SX,'TI=",F5.0,2X,'TO=’ ,F5.0,2%, PI’, F5.1,2%, P00’, 
275.1,2X,'FI=’,F4.3,2X,'FOo’, F4.3/,5X, 'ETAHC™’,F6.4,2X,'GBAR@’, 
>F6.4/,5X%,"ERPCa’ ,F6.4/) 
RETURN 
C ****SCOMPUTE LPC 
30 BIN©ENTHEN( FIN, TIN) 
PHIINSENTREN( FIN, TIN) 
PHIOIS=PHIIN+RGAS( FIN) /CJ*ALOG( PR) 
CALL PHIFN(FIN, PRIOIS, TOIS) 
HOIS*ENTHFN(FIN, TOIS) 
HO@HIN+ (HOIS-HIN) /ETALC 
WORKLC=RO-HIN 
CALL HFN(FIN,&#0,TO) 
PHIO@ENTREN( FIN, TO) 
FO-FIN 
PO=PIN¢PR 
GBAR=GAVG(FO, TIN, TOIS, PRIIN, PRIOIS) 
c WRITE(6,*)'TO=’,TO,TIN]@', TIN, 'PO=’,PO,'PIN@=', 
c >PIN, ‘ETALC#", ETALC, 'GBAR@=',GBAR 
TEMPLCL©ALOG(TO/TIN) 
TEMPLC2=AL0G(PO/PIN) 
ELPC-(GBAR-2.0)/GBAR*TEMPLC2/TEMPLC1 
¢c WRITE(6,*)ELPC=’,ELPC 
WRITE(9,40) 
WRITE(9,50)TIN,TO, PIN, PO, FIN, £0, 8TALC,GBAR.ELPC 


aaAAN 
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40 FORMAT(5X,’STATION IDENT: LP COMPRESSOR, STATIONS 1-2'/) 

50 FORMAT(5X,‘T1=",F5.0,2X,'TO#’,€5,0,2%, Pre’, F5.1,2K,'POo", 
>F5.1,2x,' Fie’, F4.3,2%,° FO’, F4.3/,5%,'RTALCH’ ,F6.4,2%,'GBAR=", 
>76.4/,5X,ELPC@’,P6.4/) 

RETURN 
END 


C tee RROPHPHRHRAAEERAAARARHEERAREREEARAEEERAEEEEEEEEERERERARORO ERED 
C tae SPAAWAVE ROTOR SUBROUTINES ARR PARR ARERRRHhENODEARAEHHRAReHeen 
C tee HRORRERERARRAREREARHEEEERHAREREEERHEREEAAOREEAEEHEARA EERE OD 


SUBROUTINE ROTOR(PRATI,PT!,PIN,PO,TIN, TO, PREF,RREF,TRE?, PIN, FO, 

>SPR,TPP,AMEX,TCCIN, IDB1, CONFIG) 
COMMON/LOSSES/DIFLOS,CCLOSS ,ABLOSS, BPLOSS 
COMMON /BURN/HFUEL, PLAV, TMAXCC, TMAXAB , THAXWR 
COMMON/CONST/CJ,G 
puMMyeO. 
GGUESS#1.3 
FOeFIN 
TSO@TMAXCC/(1.4(GGUESS-1.)/2.*AMEX**2,} 
KCT=0 

10 KCT#KCT+41 

HSO@ENTHFN(FO,TSO) 
GLOCO#GLOC(FO,TSO, ASO) 
PO@PIN*TPR*(1.-CCLOSS) 
PSO=P0/(1.4+(GLOCO-1.)/2.*AMEX**2.)#*(GLOCO/(GLOCO-1.)) 
PREF<PSO/SPR 


Con----- CONVERT DENSITY FROM kg/m*? to lba/ft*3-~----~------------~- 


TREF«PREF/RGAS( FO) /RREF*2306.6587 

STR@SPR**( (GLOCO-1.)/GLOCO) 

TSO=STR*TREF 

TO=TSO*(2.4+(GLOCO-1, ) *AMEX**2./2.) 

WRITE(‘,*)'TREFe’, TREE, 'TMAXWR=', TMAXWR, 'KCT=’, KCT, ’FO=',FO 

1F((ABS( TREF-TMAXWR) .LT.5.0).OR.(KCT.GT.5)}) GOTO 20 

TMAXWR=TREF 

re erese COMPUTE TEMP INTO THE BURNER (TCCIN}~-------~----~--------~-- 

TCCINeTIN®(PRAT1*PTI/PIN) **#((GLOCO-1.)/GLOCO) 
WRITE(6,*)'PRATIl<=’,PRATI,’ PTle',PTl 
WRITE(6,*)’TCCIN = ',TCCIN 

CALL BURNER (DUMMY, TCCIN, DUMMY, TMAXWR, FIN, FO, IDB], KODE,1) 

GOTO 10 

20 CONFIG=0 

IDBieIDBI1+1 

WRITE(6,*)°IDBl=', IDB) 

WRITE(9,30) 

WR1sze(9,40)TCCIN, TO, TREF, PIN, PO, PREF, FIN, FO,GLOCO,IDB1 

30 FORMAT(5X,'STATION IDENT: WAVE ROTOR, STATIONS 3-4'/) 

40 FORMAT(SX,’'TCCIN|',#7.1,2X,'TO*', F7.1,2X, 'TREF=’ ,FI.1/,5%, 
>'PIN]’,76.1,2X%,°PO~’,f6.1,2X,' PREP =’, F6.1/,5X,'FIN@', F6.4, 2%, 
>'FO#',F6.4,2X,GLOCO=’ ,F5.3,2X, 'IDBl=’ ,I1/) 

RETURN 
END 


aee PHARAOH EMP OHEDRARHORPADDROREREREHERARRHAARERRHERAHERAEPR OED 

SUBROUTINE BURNER( PIN, TIN, PO,TO,FIN, FO, 1DB2, KODE, CONFIG) 
INTEGER CONFIG 
COMMON /ETAS/ETALC, ETAHC, ETACC, ETALT, ETAAS,VCOEFP 
COMMON/LOSSES/DIFLOS,CCLOSS,ABLOSS,BPLOSS 
COMMON/BURN,/RFUEL, FLAV, TMAXCC, TMAXAB, TMAXWR 

pa as CHECK COUNTER OFERATION 

aa a aaa ar CONFIG=1 GOTO WAVE ROTOR CALC 

SersssoSes=sese- CONFIG RESET = 0 IN ROTOR SUBROUTINE 
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warren ween nen IDB INCREASED BY 1 IN ROTOR SUBR 


----CONFIG=0 AND IDBel GOTO MAIN CC CALC 





mer eee ewe eersern CONFIG=0 AND IDB>1 GOTO A/S CALC 


IF(CONFIG.2Q.1) GOTO 50 
IDB2=IDB2+1 
WRITE(6,*)'CONFIG=’,CONFIG,® IDB2=', IDB2 
IF(1IDB2.£0.1) GOTO 10 
GOTO 20 


meee MAIN BURNER CALC------~----------~~-~-~---------~-~-~----------- 


10 TO=TMAXCC 
ETAB*ETACC 
DELTAF=(A(TO)-A(TIN)4PFIN¢(B(TO)-B(TIN)))/(MFUEL+ETAB* FLAV-B (TO) ) 
FOeDELTAF+FIN 
PO=PIN*(1.-CCLOSS) 
HIN-ENTHFN( FIN, TIN) 
HOwENTHFN( FO, TO) 
PHIIN©ENTRFEN( FIN, TIN) 
PHIOSENTREN( FO, TO) 
WRITF(9,60) 
WRITE(9,70)TIN,TO,PIN, PO, FIN, FO 
RETURN 


----- AFTERBURNER CALC-----~--~--------------------~------~------------ 
----- A/B ON 


20 IF(KODE.£Q.0) GOTO 30 
TO=TMAZAB 
ETABeETAAB 
DELTAFe(A(TO)-A(TIN)+FIN*(B(TO)-B( TIN) ))/( HFUEL+ETAB*FLAV-B( TO) ) 
GOTO 40 


wnoee A/B OFF 


30 TOeTIN 
DELTAF=0. 
40 FO=§DELTAF+FIN 
PO=PIN*(1.-ABLOSS) 
RIN@FNTHFN( PIN, TIN) 
RO-ENTHFN(FO,TO) 
PHIIN@ENTREN( FIN, TIN) 
PRIO“ENSREN( FO, TO) 
WRITE(9,80) 
WRITE(9,70)TIN,TO,PIN, PO, FIN, FO 
RETURN 
--WAVE ROTOR CALC------ een en et er er meee ent m nme eoses coe eressereee 
50 TO=TMAXWR 
ETAB<ETACC 
DELTAF=(A(TO)-A(TIN)+FIN*(B(TO)-B( TIN) ))/( HFUEL+ETAB*FLAV-B( TO) ) 
FOeFIN+DELTAF 
60 FORMAT(5X,’STATION IDENT: MAIN BURNER, STATIONS 3-4'/) 
70 FORMAT(SX,‘TI=',F5.0,2X,'TO=', F5.0,2X,'Pi=’,P5.1,2X,'"PO='’,P5.1, 
>2K,° FI, F4.3,2K,'°FOw! 4.9 /) 
60 FOKHAT( OX, STATION IDENT: AFTERBURNER, STATIONS 9-10'/) 
RETURN 
END 


¢ ana RARRARARORARAREARARERARAEERADEEREREAREEKRERERADAEAAREDAARAEARERA 


SUBROUTINE MIXCP(FINI,FIN2,PIN1,TIN1,PIN2,TIN2,FO,PO,TO,VAR,IOM) 
COMMON/ETAS /ETALC, ETAHC, ETACC, ETAHT, ETALT, ETAAB, VCOEFF 
COMMON/LOSSES/DIFLOS,CCLOSS,ABLOS8, BPLOSS 
IDMeICN+1 
IF(TDM.£Q.1) GOTO 10 
ZETA2VAR 
GOTO 20 
10 ZETA@VAR,‘(1.-VAR) 


60 














aaa 


----CALCULATE OPTIMUM BYPASS RATIO 


























































20 HINL©ENTHEN( FINI, TINI) 
HIN2=ENTHFN( FIN2, TIN2) 
PHIINI©ENTREN(FINI,TIN1) 

PHIIN2Z@ENTREN( FIN2, TIN2) 
BO=(ZETA*(1.4FIN1) *HINL+(1.4+FIN2)*HIN2)/( ZETA*(1.¢FIN1)4(1.4PIN2)) 
FO=(ZETA*FINI+FIN2)/(1.+ZETA) 
PO=PIN2 
CALL HFN(FO,80,TO) 
PHIOwENTRFN(FO,TO) 
IF(IDM.£Q.1) GOTO 30 
GOTO 40 
30 WRITE(9,50) 
WRITE(9,60)TINI,TIN2,T0,PO,FINI,FIN2,FO 
RETURN 
40 WRITE(9,70) 
WRITE(9,60)TIN1,TIN2,TO,PO, FINI, FIN2,FO 

50 FORMAT(5X,'STATION IDENT: PRE LPT MIXING, STATIONS 3,5-6'/) 

60 FORMAT(5X, ‘Tl=',F5.0,2xK,'T2=',f5.0,2X,'TOm’,F5.0,2X,'POw’, F5.1, 

>2K,'Fle’,F4.3,2%,°F2=" ,F4.3,2%,'FOo’, F4.3/) 

70 FORMAT(SX,’STATION IDENT: PRE NOZZLE MIXING, STATIONS 9,10-11'/) 
RETURN 
END 


C tte PRARAEHRARRARRRERARARERARRERRERERARAERARAARARARHARRARAREERERHRAHAHERER 


SUBROUTINE TURBIN(PIN,TIN,PO,TO,TOIS, FIN, FO,HO,PRAT,PRLC, PRAT, PRLT 
>,X1,BETA,RASET,IDT) 
COMMON /ETAS /ETALC, ETAHC, ETACC, ETAHT, ETALT, ETAAB, VCOEFF 
COMMON/LOSSES/DIFLOS,CCLOSS,ABLOSS,BPLOSS 
COMMON /WORK1 /WORKLC, WORKHC 
COMMON/CONST/CJI,G 
IDT#1DT+1 
RIN@ENTHFN( FIN, TIN) 
PHIIN@ENTREN( PIN, TIN) 
If(IDT.£Q.1) GOTO 10 
IF(KASET.EQ.1) GOTO 20 
GOTO 30 


----HPT CALCULATION-----~-----------------~------------------- 


10 HO*HIN-WORKHC/((1.-X1)*(1.4+FIN)) 

CALL HFN(FIN,HO,TO) 

PRHIO@ENTRFN( FIN, TO) 

ROISeHIN-(HIN-HO) /ETART 

CALL BFN(FIN.#OIS,TOIS) 

PHIOIS@ENTRIN(FIN, TOIS) 

PRAT#=1./EXP(CJ/RGAS( FIN) *(PHIOIS-PHIIN)) 

FO=FIN 

PO=PIN/PRAT 

GBAR=GAVG( FO, TIN, TOIS, PHIIN, PHIOIS) 
TEMPHT1L©ALOG( PO/PIN) 
TEMPHT2=©ALOG(TO/TIN) 
EHPT=@GBAR/(GBAR-1.)*TEMPBT2,/TEMPHRTL 
WRITE(6,*%)'TO=',TO,’TIN]@’,TIN,’PO=’,PO, 

>’ PINe', PIN, ’GBAR=',GBAR, "ETAHT=',ETART 
WRITE(6,*)’EHPT=" ,fHPT 
WRITE(9,40) - 
WRITE(9,50)TIN,TO, PIN, PO, PIN, FO, ETAHT,GBAR, SHPT 

RETURN 


20 PRLT#1./((1.-BPLOSS)/(1.-CCLOSS) *PRLC/PRAT®*PRHT} 
PHIOIS@ENTRFEN( FIN, TIN) +RGAS( FIN) /CJ*ALOG(1./PRLT) 
CALL PHIFN(FIN, PRIOIS, TOIB) 

BOIG@ENTHFN( FIN, TOIB) 





c 

Canun 
30 

c 

c 

Cc 
40 
50 
60 

C tee 
10 
20 

C set 

c 


RO-RIN-ETALT* (HIN-HOIS) 
CALL AFN(FIN,#O,TO) 
PRIOsENTREN( FIN, TO) 
BETA©(1.+¢FIN)*(HIN-RO) /(WORKLC)-1. 
WRITE(6,*)*TURBIN BETA=' ,BETA 
FO=FIN 
PO-sPIN/PRLT 
GBAR-GAVG( FO, TIN, TOIS, PHIIN, PRIOIS) 
TEMPLT1-ALOG( PO/PIN) 
TEMPLT2-ALOG(TO/TIN} 
ELFT=GBAR/(GBAR-1.)*TEMPLT2/TEMPLT1 
WRITE(9,60) 
WAITE(9,50)TIN, TO, PIN, PO, FIN, FO, ETALT,GBAR, ELPT 
RETURN 
-LPT CALCULATION BYPASS PRESCRIBED-~---~-~--------~-------- 
HO=jHIN-(1.4+BETA) *(WORKLC)/(1.4FIN) 
CALL HFN(FIN,HO,TO) 
PELO=ENTREN( FIN, TO) 
HOIS@HIN+(HO-HIN) /ETALT 
CALL RFN(FIN,HOIS,TOIS) 
PRIOIS*ENTREN( FIN, TOIS) 
PRLT=1./EXP(CJ/RGAS(FIN)* (PHIOIS-PHIIN) ) 
TO-rINn 
PO=PIN/PRLT 
GBAP~GAVG(FO, TIN, TOIS, PHIIN, PHIOIS) 
TEMPLTL=ALOG( PO/PIN) 
TEMPLT2=<ALOG(TO/TIN) 
F.UPT=GBAR/(GBAR-1.) *TEMPLT2/TEMPLTL 
WRITE(6,*)*TO=’,TO,*TIN]=’, TIN, 'PO=’,PO,’PIN=’,PIN, 
>'GBAR=’,GBAR, ETALTe’ , £ETALT 
WRITE(6,*)ELPT=',ELPT 
WRITE(9,60) 
WEITE(9,50)TIN, TO, PIN, PO, FIN, FO,ETALT,GBAR, ELPT 
FORMAT(SX,’STATION IDENT: BP TURBINE, STATIONS 4-5'/) 
FORMAT(SX,*T1=©',F5.0,2X,°TO=’, F5.0,2K,Pl=’,F5.1,2X, POW’, 
>P5.1.2z,°FI—’, F4.9,2K,' FO’, F4.3/,5X, ETA=', F6.4,2X, 'GBAR=", 
>f6.47,5%,'UT@',F6.4/) 
FORMAT(SX,'STATION IDENT: LP TURBINE, STATIONS 6-8'/) 
RETURN 
END 
ARERAERARHAARRARERARERERARARREERAAAAARRAREHRHRARARARAERAERERRAREA 
SUBROUTINE DUCT(PIN, TIN, FIN,PO,TO, FO, RO) 
COMMON/LOSSES/DIFLOS,CCLOSS, ABLOSS, BPLOSS 
PO-PIN*(1.-BPLOSS) 
TO=TIN 
ro=-0.0 
RIN PNTHFN( FIN, TIN) 
PULIR|OZNTPEN( TIN, TIN) 
HO@ENTHFN(FO,TO) 
PRIO@ENTRFEN( FO, TO) 
WRITE(9,10} 
WRITE(9,20)TIN,TO,PIN,PO,FIN,fO 
FORMAT(SX,*STATION IDENT: BYPASS DUCT, STATIONS 2-7'/) 
PORMAT(S¥, °TI=’,F5.0,2K,'TO~’ ,f5.0,2K, Pim’ ,F5.1,2K,'PO=’, 
>F5.1,2x,°’Fr]’,F4.3,2X,'FO3=’,F4.3/) 
RETURN 
END 
RERPHRAARARAARARARARERERERARARAHRARAAEARERHRERARARHEARARARAREARRADEROR 


BUBROUTINE MIXCA(PIN2, PIN), TIN2, TIN] ,BIN2,HINI,FIN2, INL, BETA,CMIN 
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>1,PIMAX,TO,FO,PO,CMO) 

COMMON/CONST/CJ,G 
CHRAAARAARAHORARREHAREARERERERREEREEESEAARESTATEMPENT FUNCTION DEFINITIONS 

PHI FN({ GAMMA, AMACH) @AMACH*#2,®(1.4(GAMMA-1.)/2.*AMACH**2.)/((1.4GA 

>MMA* AMACH*#2,)*#2,) 

PSILFN(G,GAMMA,R,AMACH)=(AMACH*SQRT(G*GAMMA/R) )/(1.4(GAMMA-1.)/2.¢ 

>AMACH**2,)**( (GAMMA41.)/(2.*(GAMMA-1.))) 
CAERMRADHHIRAAAAAREARRAERRARRHAEHARRARAAREBOLYVE FOR FAN STREAM MACH NUMBER 


c CHECK IF BYPASS GREATER THAN 0.0 
IF (BETA.GT.0.0) THEN 
GOTO 5 

ELSE 
TO=TINI 
PO-PINI *PIMAX 
RO-HINI 
FO=FIN1 
CMO-CMINI 
AI2A11=0.0 

ENDIF 

Go TO 9 


CAAHRAHRAARERAARAAAARARARAHERRORARARRAAD 
5 GLOCI 2=GLOC( FIN2, TIN2, RIN2) 
c WRITE(6,*) AMACHB~’ ,CMINI 
GLOCI1@GLOC( FINI, TINI,RIN1) 
TLAGL=(PIN2Z/PIN1*(1.4(GLOCI1-1.)/2.*CMINL®*2.)#*#(GLOCI1/(GLOCI 
>1-1.)))**((GLOCI2-1.)/GLOCI2) 
Cc WRITE(6,*)‘ FLAG1@=", FLAG) 
IF (BETA.GT.0.0) THEN 
If(FLAG!.LE.1.0) GOTO 60 
ENCIF 
CMIN2=SQRT(2./(GLOCI2-1. )#((PIN2Z/PINI*(1.4(GLOCII-1.)/2.*CMINI**2. 
>)**(GLOCIL/(GLOCII-1.)))*4( (GLOCI2-1.)/SLOCI2)-1.)) 
IFf(CMIN2.GE.0.8) GOTO 40 
CAAAHAKAHHAAARRAHEREHHRERAERAEHRAESOLVE FOR MIXER INLET AREA RATIO 
AI2A11L@BETA/(1.4FIN1) *CMINI/CMIN2*SQRT( (GLOCI1I*RGAS(FIN2)*TIN2Z¢(1. 
>+(GLOCI1-1.)/2.*CMIN1**2.))/(GLOCI2*RGAS(FIN1)*TIN1*(1.4(GLOCI2-1. 
>)/2.*CMINZ9*2.))) 
CHAAHAHAHRHAREAEAREREAAAHERAAHERERHREREOLVE FOR MIXER TEMPERATURE RATIO 
RIN2@RGAS(FIN2) 
RIN1=RGAS(FIN1) 
CPIN2=CP( FIN2, TINZ) 
CPINI=CP(FINI,TINI) 
BETAPR=BETA/(1.4+FIN1) 
CPO=(CPIN14+BETAPR*CPIN2)/(1.+BETAPR) 
TOeTIN1*(CPINI/CPO*(1.+BETAPR*CPIN2*TIN2/(CPINI*TINI))/(1.+BETAPR) 
>) 
TAUMIX=TO/TINI 
CHRARARHAAAEDARHOERREARAHARAPAARARREQOLVE FOR MIXER EXIT MACH NUMBER 
FO= (BETAS IN2+FIN1)/(1.4BETA) 
HO-ENTHFN( FO, TO) 
RO=RGAS( FO) 
GLOCO=GLOC( FO, TO, HO) 
PHIINZ=PHIJFN(GLOCI2,CMIN2) °- 
PHIINI@PHIJEN(GLOCI3,CMIN}) 
PHIOe(((1.4+BETAPR)/(SQRT(1./PHIIN1)+BETAPR*SQRT(1./PHIIN2) *SQRT(RI 
2N2*TIN2Z*GLOCIL/{ RINL*TINI*GLOCI2))))**2.)*GLOCI1*RO* TAUMIX/(GLOCO* 
>RINL) 
CmO*SQRT(2.*PHIO/({(1.-2.*GLOCO*PHIO)+SORT(1.-2.*(GLOCO41.)*PHIO))) 
CARRREAEHORHAHERERADEDEREHHERO RADAR HRESHHEOLVE FOR MIXER PRESSURE RATIO 
POI@PINI*(1.4BETAPR)*(1./(1.4AI2AI1)) *SQRT( TAUMIX) *PSILFN(G,GLOCIL 
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c 
c 


aaAN 





>, RIN1,CMINIL)/PSI1FN(G,GLOCO,RO, CHO) 
A1@PSIIFN(G,GLOCI1,RIN1,CMINI) 
A2@PSI1FN(G,GLOCO,RO, CMO) 
PIMIXI©POI/PINI 
PIMIX=PIMIXI*PIMAX 
PO=PIMIX*PINI 
9 WRITE(9,10) 
WRITE(9,20)TIN1,TIN2,TO, PIN, PIN2,PO 
WRITE(9,30)FIN1,PFIN2,FO,CMINL,CMIN2,CMO,AI2ZAI1 
RETURN 
10 FORMAT(SX,’STATION IDENT: CONST. AREA MIXER, STATIONS 7,86-9'/) 
20 FORMAT(5X,‘'T1=',F5.0,2X,°T2=",F5.0,2K,'TOm’,F5.0,2%, Plo’ ,P5.1,2%, 
>*Plm’,FS5.1,2X,"°PO"’,F5.1/) 
30 FORMAT(5X, lao’, 04.3,2X,'F2=', F4.35,2%, 'FOm’, F4.3,2K, Mla’, P4.2,2K, 
> Mle’, F4.2,2K,'MOo’, F4.2,2K%, ARAT@" ,F4.2/) 


40 WRITE(6,50) 
50 FORMAT(5X,’FAN STREAM MACH NO IS GREATER THAN/EQUAL TO 0.6’) 
RETURN 


60 WRITE(6,70) 
70 FORMAT(5X,'FAN STREAM MACK NO I8 LESS THAN/EQUAL TO 0.0’) 
RETURN 
END 
REA RARAAERAROREEEEAERERAREEAEREEEARRERARARAERARERARDHADRERORRERDR EDEL 


BUBROUTINE NOZZLE( PIN, TIN, PO,TO,TOIS, FIN, PRLC, VEXIT, EXMACH, PSO, 
>IFLAG) 
COMMON/ETAS/ETALC, ETAHC, ETALC, ETAHT, ETALT, ETAAB, VCOEFF 
COMMON/LOSSES/DIFLOS, CCLOSS,ABLOSS, BPLOSS 
COMMON/CONST/CJ,G 
PRN@1./((1.-DIFLOS)*#(1.-BPLOSS)*(1.-ABLOSS) *( PRLC) ) 
PRN@PSO/PIN 
HIN@ENTHFN( FIN, TIN) 
PHIIN@ENTREN( FIN, TIN) 
GLOCO“GLOC( FIN, TIN, HIN) 
PRN=PSO/PIN 
IF(PRN.GE.1.0) GOTO 30 
write(6,*) "GLOCO~’ ,GLOCO 
write(6,*)'PSO='’,PSO 
weite(6,%)'PRN=’,PRN 
PRIOIS~PHIIN+RGAS( FIN) /CJ*ALOG( PRN} 
CALL PHIFN( FIN, PHIOIS, TOIS) 
HOIS-ENTHFN( FIN, TOIS) 
PO=PIN*PRN 
ETAN=VCOEFF##2, 
HO=HIN-ETAN*(HIN-HOIS) 
CALL HFN(FIN,HO,TO) 
VEXIT=SQRT(2.*G*CJ*( HIN~HO) } 
AEXIT=SQRT(G*RGAS( FIN) *TO/(1.-RGAS( FIN) /CJ*TO/HO) } 
EBXMACH@VEXIT/AEKIT 
WRITE(9,10) 
WRITE(9,20)TIN, TO, PIN, PO, PIN, VEXIT, EXMACH 
WRITE(9,25)GLOCO, PRN, AEXIT 
GOTO 40 
10 FORMAT(5X,"STATION IDENT: EXHAUST NOZZLE, STATIONS 11~-12°/) 
20 FORMAT(5X,’TI©’,F5.0,2X,’TO=',F5.0,2%, Pleo’, F5.1,2%,'PO~’, F5.1, 
>2X, °F Ie’, F4,3,2K,'VO0’, F5.0,2%,'MO=', F4.2) 
25 FORMAT(5X,"GLOCO=’,f6.4,2X,'PRN=’,F6.5,2X, "AEXIT=",¥6.1/) 
320 IFLAG-1 
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40 RETURN 
END 
¢ RHR HRHEAARARAARARARARARARARARARRARARRAERAEHRHERERAAARHERRERRARARAREREREREREE 
c 
> SUBROUTINE PERF(AMACHS, FIN, VO, VEXIT, BETA, SPTHR, SFC, EXMACR) 
COMMON/CONST/CJ,G 
SPTHR|((2,4 FIN) *VEXIT-VO)/G 
SFC=FIN*3600./SPTHR 
WRITE(9,10) 
WRITE(9,20)AMACH8, SPTHR, SFC, BETA, EXMACH 
WRITE(9,25) 
10 FORMAT(5X,* PERFORMANCE PARAMETERS’ /) 
20 FORMAT(5X, "AMACHS=',¥6.3/,5X,'SP. THRUST =’, f6.2,2X,°SFC = °,F5.3 
>,4X%,'BETA= °,F6.3,2X, 'EXMACH@=',F6.3) 
25 FORMAT(5X,°(LBF/LBM/S)',6X,'(LBM/HR/LBF )'/) 
RETURN 
END 
Cc RAR HRARHARHRARHRARARAAARERARRAARARANARARARHARARARARARRARARAHARRAHREARHERKRERKRE 
c 
SUBROUTINE HFN( FRATIO, HEX, TEMPEX) 
DUMMYel, 
ITERNO“0 
IHEXeIFIX(HEX/100.)*100 
HONE@ FLOAT ( LHEX) 
HTWO<HONE+100. 
CALL TEMP( HONE, ATWO, DUMMY, DUMMY, TEMP1, TEMP 2) 
10 ITERNO-“ITERNO+1 
IF(ITERNO.GT.50) GOTO 20 
HONEENTHFN( FRATIC, TEMP1) 
HTWOeENTHFN( FRATIO, TEMP2} 
SLOPE~( HTWO-HONE) /( TEMP2-TEMP1) 
TEMP1=TEMP1+(HEX-HONE ) /SLOPE 
HONZ@ENTHFN( FRATIO, TEMP) ) 
If (ABS (HONE-HEX).GT.0.002) GOTO 10 
TEMPEX=-TEMP1 
- GOTO 30 
20 WRITF‘6,40) 
40 FORMAT(5X,'SECANT METHOD ITERATION DOES NOT CONVERGE’ /} 
30 RETURN 
END 
Cc RAR RARARAARAARARARRERARAERHARRRERARARARRRARARAARRARADARHARARARHARARAEREE 


SUBROUTINE PHIFN( FRATIO, PHIEX, TEMPEX) 
DUMMY=1, 
ITERNO=0 
IPHIEX@IFIXN( PHIEX*10. ) 
PHIONE=FLOAT( IPRIEX)*0.1 
PHITWO-<PHIONE+0.1 
CALL TEMP{ DUMMY, DUMMY, PHIONE, PHITWO, TEMP], TEMP2) 
10 ITERNO“ITERNO+1 
IF(ITERNO.GT.55) GOTO 20 
PHIONE<ENTRFN( FRATIO, TEMP) } 
PHITWOSENTRFN( FRATIO, TEMP 2} 
SLOPE=(PHITWO-PHIONE) /( TEMP2-TEMP1) 
TEMPL=TEMP14( PHIEX-PRYONE) /SLOPE 
PHIONE@ENTRI?2( FRATIO, TEMPL) 
If (ABS ( PHIONE-PRIEX).GT.0.000003) GoTO 10 
TEMPEX=TEMPL 
GoTo 40 
20 WRITE(6,35) 
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10 
20 


FORMAT(5X,’BSECANT METHOD ITERATION DOES NOT CONVERGE’ /) 
RETURN 
END 
Pee SEU VUES PSEC IPS IriVirriririiiritiiirrriitri tt itil isrii iti i ty 
FUNCTION STATEMENTS 
RORADREARAREEAERERARREA DARE ORAERAREREREREREEEHAERARERERERARERERE 
FUNCTION ENTHFN( FRATIOC, TEMP) 
IP(TEMP.LT.0.0) GOTO 10 
€120.24062 
C220.0177248-03 
C3=0.0380562-06 
C420 .0126622-09 
C5=0.00130122-12 
D1=0.22091 
D2-0,51822£-03 
D3=“0.194626-06 
D4-0.0450898-09 
D5=-0,00432758-12 
ENTHFNe(1./(1.+F RATIO) )*(CI*TEMP-C2/2.*TEMP**#2 4037/3. 9° TEMPO*? ~Ca/ 
>4.*TEMP**4.4+C5/5.*TEMP**5. )+(FRATIO/(1.4¢F RATIO) )*(D1*TEMP4+D2/2.°TE 
PMPO*2.-03/3.*TEMP4**3 .4D4/4. *TEMP**4,-D5/5.*TEMP**5.) 
GOTO 25 
WRITE(6,20)TEMP 
FORMAT(SX, ‘NEGATIVE TEMP IN ENTHFH ENCOUNTERED=’,F6.1/} 
CONTINUE 
RETURN 
END 


FUNCTION ENTREN( FRATIO, TEMP) 
IF(TEMP.LE.0.0) GOTO 10 
C1-0.24062 
C2-0.0177248-03 
C3-0.038056E-06 
C4=0.012662£-09 
c5=0.00130128-12 
Die-0.22091 
D2=-0.51022£-03 
D3I-0.19462F-06 
D4-0.045089E-09 
05-0.00432758-12 
ENTREN@=(1./(1.4F RATIO) }#(CL*ALOG( TEMP )-C2*TEMP+C3*0. S*TEMP#*#2,-C4/ 
>U.*TEMP**3.4C5/4.*TEMP**4, )4( FRATIO/(1.4+FRATIO)) *(DI*ALOG( TEMP) 4+D2 
>*TEMP-D3*0.5*TEMP**2.4D4/3.*TEMP**3.-DS/4.*TEMPOed.) 
RETURN 
WRITE(6,20) TEMP 
FORMAT(5X,NEGATIVE/ZERO TEMP IN ENTREN ENCOUNTERED =',F6.1) 
RETURN 
END 


FUNCTION A(T) 
Cl=,24062 
Cl201.7724£-5 
C33. 8056£-6 . 
Cae). 2662E-11 
CS—1.3012E-15 
dle. 22091 
D2-5.1822e-4 
D3-1.9462E-7 
p4-4,50898-11 
D5=4.32752-15 
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AwC1L*T-C2/2.9T#2 4039/3. OTH83 CAG. OTHHE 4055/5 HT ANS, 
RETURN 
END 


FUNCTION B(T) 
€i=.24062 
C2~1.7724£-5§ 
€3-3.8056E-8 
C4=1.2662£-11 
CS=1.3012£-15 
Dl~.22091 
D2—8.18Z2e-4 
D3=1.9462E-7 
D4=4.5069E-11 
DSee.3275E-15 
BeDL*T4D2/2.8TH#*2 .-D3/3.*THh3 404 /4, eTH94.-D5/5.eTH*8, 
RETURN 

END 


FUNCTION CP(FRATIO, TEMP) 

C190,24062 

C2=0.017724E-03 

C320 CIAOSEE-06 

C4=0.912662E-09 

C5-0.0013012z-12 

Diet. 22091 

D290 .£1822£-03 

D3=-0.194628-06 

D4=7.9459898-09 

D5=0,0943275E-12 

CPa(L./(1.4+FRATIO) )*(C1-~CZ*TEMP+C3*TEMP**2,-CA*TEMP**3 .4CS*TEMP AHA 
Spr CER DTAO/ LL so PRRTAO VET (OLS OATS ERE PAST ERENT a: SDASTERE Seo <cDer zen 
persed.) 

RE AUN 

END 


FUNCTION RGAS( FRAT) 

WMOL + (1 + FRAT! /(0.03452240.035648*FRAT) 
RGIS 1545.43 /WMOL 

RETURN 

END 


FUNCTION GAVG(FRAT, TIN, TOIS, PHIIN, PHIOIS} 

COMMON/CONST/CI.G 

WRITE(6,5) FRAT, TIN, TOIS, PRIIN, PHIOIS,CJ 
FORMAT(5X,F3.1,5X,2F8.1,5%,2F7.2,5X,F6.2} 
GAVGe1./(1.~(RGAS( FRAT) /CJ)*(ALOG(TOIS/TIN) /(PHIOIS-PHIIN))} 
PrTeN 

ériv 


FUNCTION GLOC(FLOC,TLOC,HLOC) 
COMMON/CONST/CJI,G 
GLOTe1.7(1.-RGAS(FLOC)/CJ*TLOC/BLOC) 
RETURN 

END 


SUBROUTIN® TEMP(HONE, TWO, PRIONE, PHITWO, TENP], TEMP2) 
JF((HONE.%Q.0.00).OR.(PHIONE.LE.1.3)) GOTO 10 
IF((RONE.EQ.200.).OR. (PHIONE.£ZQ.1.4)) GOTO 20 
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IF ((RONE.2Q.200.).OR. (PHIONE.£Q.1.5)) GOTO 30 
If ( (HONE.2Q.300.).OR.(PHIONE.EQ.1.6)) GOTO 40 
If ( (HONE. £Q.400.).OR.( PRIONE.£Q.1.7}) GOTO 50 
IP ((HONE.£EQ.500.).OR.(PHIONE.£Q.1.8)) GOTO 60 
IF( (HONE.GE.600.).OR.(PRIONE.GE.1.9)) GOTO 70 


WRITE(6,1)HONE, PHIONE 


FORMAT(5X, NOTHING MATCHED: HONE@' ,F13.9,5X,*PRIONE=" ,F13.9/) 
STOP 
TEMPL«0, 
TEMP2=450. 
RETURN 
TEMP1=300. 
TEMP2=800. 
RETURN 
TEMP1=S500. 
TEMP2=1250. 
RETURN 
TEMP1=750. 
TEMP2=1600. 
RETURN 
TEMP1#1100. 
TEMP2=2000. 
RETURN 
TEMP1-1500, 
TEMP2=2400. 
RETURN 
TEMP1=2200. 
TEMF2~3300. 
RETURN 
END 
RAARAAARAARADAHARHRARAERARAAARAARARAHERARAARAAREEHHEARHRAREERARERRERARAD 
PLOT SUBROUTINES USING "DISSPLA’ SOFTWARE - FACILITY DEPENDENT 
ARARHHARRARAARARRAREHERERRARARAREARAAAAEHRARAHARARHRARARAAERERRKRERA 
SUBROUTINE PLOT( XARAY, YARAY, ZARAY,BPR,JJ,11) 
DIMENSION XARAY(1I1),YARAY(I1),ZARAY(I1),BPR(11) 
PEAL XARAY, YARAY, ZARAY,BPR 
INTEGER 11,90 
CALL PHYSOR(0.5,0.5) 
CALL AREA2D(10.,7.5) 
CALL FRAME 
CALL GRAF(0.,’SCALE’,140.,0.,"’SCALE’,100.) 
CALL GRAF(50.,'SCALE’,120.,0.6,’ SCALE’ ,1.6) 
CALL XNAME('SP.THRUST’ ,9) 
CALL YNAME(’TSFC’,4) 
GOTO(10,20,30,40,50,10,20,30,40,50),05 
CALL SETCLR(’ YELLOW’ ) 
CALL MARKER(15) 
GOTO 60 
CALL SETCLR{ ‘CYAN’ ) 
CALL MARKER( 18) 
GOTO 60 
CALL SETCLi( MAGENTA’ ) 
CALL MARKER(5) * 
GOTO 60 
CALi SETCLR( GREEN’ ) 
CALL MARRER(9) 
GOTO 60 
CALL SETCLR{’RED’) 
CALL MARRER( 2) 
CALL CURVE( XARAY, YARAY,11,41) 


68 





c CALL ENDGR(0) 

c RETURN 

c END 

re ree 

c PLOT ROUTINE USING GRAFRIT SOFTWARE 

Cee e268 8 8 ee C6 8886 8888 SES ESSE EEE OSHS OSS SD OE NEE SRSABS Se SeVeHsASeeen 
SUBROUTINE GRSPLOTI(X1A,Y1A,Y2A,IAN4,1CP1) 

Cowem- PLOT OF SPECIFIC THRUST AND SFC VS INPUT PARAMETER 

Co----- CHOICES OF INPUT PARAMETERS ARE: 

CHsee= 1. BYPASS RATIO 

Co<ee- 2. FAN PRESSURE RATIO 

C--5—= 3. COMPRESSOR PRESSURE RATIO 

C---~9- 4. LPT EXIT MACH NUMBER 

c 


REAL X1A(50),YIA(50),Y2A(50) 
CHARACTER KILABEL*6, YILABEL*33, Y2LABEL* 33 
CHARACTER F1IG1*7,NUMB1*3, TOPLBL*10 
Cam wennen steerer nernesmennanvarsussencecannsasuceseenmesese 


C----- SPECIFY X,¥ LEFT AND Y¥ RIGHT LABELS--~--~---- wecsweere- 


YiLASEL='A: Specific Thrust lbf/(lbm/e)’ 
Y2LABEL@"Bt SFC (lbm fvel/hr)/ibf thrust’ 
C----- CHOOSE THE X PARAMETER TO PLOT 
IF (JAN4.8£0.1) THEN 
KILABEL©’ BETA * 
ELSFIF (IAN4.2Q.2) THEN 
RILABEL='’ PRLC ’ 
ELSFIF (IAN4.£Q.3) THEN 
XILABEL@=’ PRHC ’ 
ELSEIF (IAN4.£0.4) THEN 
X1LABEL@’AMACHB’ 
ENDIF 


WRITE(6,*)' INPUT: & 
WRLITE(6,*)° 0 - AUTOMATICALLY DETERMINE XIMIN AND X1MAX' 
WR1I12,6,*)}° 1 - INPUT XIMIN ANO X1MAX’ 
READ(S,*)IANSS 
Tre ‘TANS5.EQ.1) THEN 
WRITE(6,*)" INPUT XIMIN @=3’ 
PEAD(S,*)XIMIN 
WRITE(6,*)" INPUT XIMAX e@:! 
READ(5,*)X1MAX 
ELSE 
C----- AUTCMATICALLY DETERMINE THE MIN AND MAX X VALUE 
RIMINe@XIA(1L) 
XiMAKeX1A(1) 
DO Je1l,ICPI 
RITEMPRXIA(I) 
X2TEMP@XIA( I) 
If (FIMIN.GT.XITEMP) THEN 
XIMIN@XITEMP 
ENOIF 
IF (X1MAX.LT.X2TEMP) THEN 
XIMAX@X2TEMP 
ENDIF 
exD DO 
ENDIF 


WRITE(6,*)'RIMIN]’ ,XIMIN, 'XIMAX=" , XIMAX 


aAaAaQ 
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Cwe--- INPUT THE NUMBER OF THE FIGURE TO PLOT (FIG 1)--- 
WRITE(6,*)}*‘ INPUT THE FIGURE NUMBER (fe 1 OR 2...)° 
READ(5,°(A9)* )NUMBL 
¢c WRITE(6,*)'NUMB1=' ,NUMBL 
FIGle' FIGURE ’ 
TOPLBL=FIG1//NUMB1 
Cc WRITE(6,'(1X,A11}°)’TOP LABEL =’, TOPLBL 
Conse eee ee eS ess Be SSS eS SS SSS ee essere ewaseesseeaserness 
eee CHECK THE DATA INPUT 
DO 1l#1,1CP1 
WRITE(6,*)12,X1A(11),-YIA(I1L),¥2A(11) 


owece CONSTRUCT THE PLOT 


¢ 
c 
c 
c END DO 
c 
Cc 


WRITE(6,*)*STANDBY FOR PLOT’ 


C----- OPEN GKS, OPEN AND ACTIVATE WORKSTATION. 
CALL OPNGKS(1,2,1} 

c 

c SET UP THE COLOR TABLE 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


GSCR(1,0,0.00,0.00,0.50) 
GSCR(1,1,1.00,1.00,0.00) 
GSCR(1,2,0.00,1.00,0.50) 
GSCR(1,3,0.00,1.00,1.00) 
GSCR(1,4,1.00,0.00,0.00) 
GSCR(1,5,1.00,0.00,1.00) 
GSCR(1,6,1.00,0.00,1.60) 
GSCR(1,7,1.00,1.00,1.00) 
GSCR(1,8,0.00,0.00,0.00) 


C----- SET VIEWPOINT AND DEFINE TRE RANGE OF USER'S DATA TO BE USED 


CALL 
CALL 
CALL 
CALL 


AGSETF ('GRID/LEFT.', 
AGSETF (’GRID/RIGHT.', 
AGSETF (‘GRID/TOP.’, 
AGSETF (‘GRID/BOTTOM.’ 
FOR FIRST PLOT 


-85) 
-85) 
» 15) 


C----- DEFINE THE TOP, BOTTOM, LEFT AND RIGHT LABELS 
C----- PUT BLANKS IN FIRST CURVE LABELS 


CALL 
CALL 
CALL 


CALL 
CALL 
CALL 


CALL 
CALL 
CALL 


CALL 
CALL 
CALL 


AGSETC ('LABEL/NAME.’, 
AGSETI ('LINE/NUMBER.’ 
AGSETC ('LINE/TEXT.’, 


AGSETC ('LABEL/NAME.’, 
AGSETI (°LINE/NUMBER.’ 
AGSETC ('LINE/TEXT.', 


AGSETC (’LABEL/NAME.’, 
AGSETI (‘LINE/NUMBER.’ 
AGSETC ('LINE/TEXT.', 


AGSETC (’LABEL/NAME.’, 
AGSETI (’LINE/NUMBER.’ 
AGSETC (*LINE/TEXT.', 


C---~- TURN OFF RIGHT AXIS 
CALL AGSETF (’AXIS/RIGHT/CONTROL.',0.0) 
Ce--= SET YMIN AND YMAX FOR FIRST P 
CALL AGSETF (‘°X/MINIMUM.’, XIMIN) 
CALL AGSETF ('X/MAXIMUM.’, XIMAX) 


CALL 
CALL 


AGSETF ('Y/MINIMUM.’, 
AGSETF (°Y/MAXIMUM.’, 
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'T') 
» 100) 


, ) 


"B’) 
27100) 
ns) | 


*L') 
+ 100) 


RR’) 
»-100) 
ery 


LOT 


0.0) 
120.0) 


CALL AGSTUP (X1A,1,1,I1CP1,1,Y1A,1,1,1CP1,1) 
CALL AGBACK 
CALL GSPLCI (7) 
CALL DASHDC ('S$SS$Ssssss’7 errr $$ssssssss’ 3,1) 
C------ PLOT CURVE 1: X1A VS Y1A 
CALL AGCURV (X1A,1,Y1A,1,I1CP1,-1) 


c 
Cessavenascescaaessoaezeas 
C----- SET UP FOR SECOND PLOT 
C----- DEFINE THE TOP, BOTTOM, LEFT AND RIGHT LABELS 
C----- THE LARGER LINE NUMBER WRITES ABOVE PREVIOUS LABELS 
CALL AGSETC (‘LABEL/NAME.’, ‘T’) 
CALL AGSETI (’LINE/NUMBER.’, 100) 
CALL AGSETC ('LINE/TEXT.’, TOPLBL) 
Cc 
CALL AGSETC ('LABEL/NAME.’, 'B’) 
CALL AGSETY (*LINE/NUMBER.’,-110) 
CALL AGSETC (’LINE/TEXT.',XLLABEL) 
c 


CALL AGSETC ('LABEL/NAME.’, 'L’} 
CALL AGSETI (’LINE/NUMBER.', 110) 
CALL AGSETC ('LINE/TEXT.’,YILABEL) 


CALL AGSETC (*’LABEL/NAME.’, '‘R') 
CALL AGSETI {’LINE/NUMBER.’,-100} 
CALL AGSETC ('LINE/VEA7.', YZLABZL) 

Cc TURN OFF THE TOP, BOTTOM AND LEFT AXES 
CALL AGSETF ('AXIS/TOP/CONTROL.’,0.0) 
CALL AGSETF ('AXIS/BOTTOM/CONTROL.’,0.0) 
CALL AGSETF ('AXIS/LEFT/CONTROL.',0.0) 


C--~-- TURN ON THE RIGHT AXIS 
CALL AGSETF ('AXIS/RIGHT/CONTROL.’,1.0) 
Cow--- CHANGE THE YMIN AND YMAX 


CALL AGSETF (°Y/MINIMUM.’, 0.0) 
CALL AGSETF ('Y/MAXIMUM.'’, 1.5) 
C----- TURN ON THE RIGHT AXIS NUMERICS 
CALL AGSETF ('RIGHT/NUMERIC/TYPE.' ,1.E6) 


CALL AGSTUP (X1A,1,1,1CP1,1,¥2A,1,1,1CP1,1) 
CALL AGBACK 
CALL GSPLCI (4) 
CALL DASHDC (’$$sFC$’'''',B,1) 
C--~---PLOT CURVE 2: X1A VS Y¥2A 
CALL AGCURV (X1A,1,Y2A,1,ICP1,-2) 


Car saee CALL FRAME TO ADVANCE THE FRAME 
C-----CLOSE GRKS, CLOS£ AND DEACTIVATE WORKSTATION 
CALL FRAME 


CALL CLSGKS 


Connenemncwesnsensesuene 


END 
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APPENDIX B 


PROGRAM WRCOMP 


A. PROCEDURES 
1. WRCOMP 
a. Log on to the VAX 
b. At the DCL prompt type: 


Srun WRCOMP1A........-.- runs WRCOMP program 
otonie? di faite wares select options from screen 


2. INTERFACING WITH ENGINE 
a. Run ENGINE program in either configuration. 


b. $type ENGINEI.DAT or 
Sprint ENGINE1.DAT/que=1la210 1 


c. From the ENGINE1.DAT determine: 
1) Pressure out of HPC 
2) Temperature out of the HPC 
3) Gamma at the HPC 
4) Fuel-air ratio (f) from the burner 


ad. Run WRCOMPI1A using these values as inputs and make 
a guess for RREF. 


e. If min doesn't equal mout from WRCOMP, then change 
the guess for RREF and iterate until they are 
equal. 


f. Run ENGINE1A in the wave rotor configuration with 
the correct value for RREF. 
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B. SAMPLE RESULTS 


FILE WRPERF1.OUT FROM WRCOMP1A.FOR PROGRAM: 


@ OF TIME STEPS Ke 2000 
CFLNUM=  0.6000000 

# OF RIEMANN ITER. (QSTOP)= 20 
GRID CELL WIDTH DX= 0.01 

PTO = 79.20000 psi 


PTOTIN 546064.8 Pa 

RTIN - 0.1956737 Ibm/ft*3 
RTOTIN = 3.134375 kg/m" 3 
PSEXIT1 ~ 75.54148 psi 
PSEXIT = 520840.2 Pa 

PREF] = 188.8537 (psi) 
PREF = 1392101. (Pa) 

xref = 0.1800000 B 

xcefl = 0.5965512 ft 


TIMING AND MASS FLOW BALANCE: 
INLET PORT OPFNS AT: 


Ne- 741 TTO7AL# 1.2576669 TIME= 0.0003768 
EXHAUST PORT CLOSES AT: 

Nei0l4 TTOTAL®= 1.7417805 TIME= 0.0005218 
INLET PORT CLOSES AT: 

N=1681 TTCIAL@ 3.1865280 TIME= 0.0009547 
P(195)=0.60590% P(201)=0.684240 P(203)=0.684240 
Min= 0.561596 Mout= 0.626005 SWL= 1 SWRe 1 


WAVE ROTOR PERFORMANCE MAP ANALYSIS: 
SPR= 0.40 TRH = 1.60 TTR =1.1356 FIN =0.0165 AMEXIT=0.9000 
PTO~ 79.20 TOTIN= 1056.00 GAMMA#1.371¢ RTINwO.1957 RTOTIN|3.1344 
PTE~126.72 FE - 75.54 PREF #186.65 TIL = 1964.0 TE = 1707.4 
TREF~ 2187.9 RGUESS=0,22520 RREFe 3.6073 


COMFUTER RUN TIME = 635.5391 gecé. 
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C. WRCOMP PROGRAM LISTING 

Sections which are commented out with a "C" contain 
features that have been disabled (eg. DISSPLA) or they contain 
amplifying remarks. ‘ 


Cc PROGRAM WRCOMP WITH VAN DER CORPUT SAMPLING AND SINGLE TIME STEP 
INTEGER QPRINT, QSTOP,5WL,5SWR - 
REAL MASSIN,MASOUT 
DIMENSION WNORM(12),IDIGT(12),XARRAY(100) 

DIMENSION P(203),R(203),U(203),A(203),5(203),xX(203) 

COMMON /SUBS/P,R,U,A,S,X 
COMMON/GLIMM1L/PGLIM,RGLIM,UGLIM, PL,RL,UL, PR,RR,UR,AL,AR,GL,GR-EPS 
COMMON/GLIMM2/DT,DX,XI 

COMMON/FUN1I/G, PA, RA,UA,RB,RMU 

COMMON/SAMPLE/WNORM, IDIGT 

COMMON XARRAY,N1 


c CALL COMPRS 

Cc CALL BLOWUP(1.5) 

c CALL PAGE(11.0,8.5) 

c CALL HWSCAL(‘'SCREEN’ ) 

Cc DATA K,SWL,SWR/2000,1,2/ 

c DATA N,CFLNUM, TTOTAL/0,0.60,0.0/ 

C DATA PTOTIN, RTOTIN, PSEXIT/547292.900,3.250,497390.1/ 

Cc DATA PREF,RREF,XREF/1243475.4,3.60,0.1800/ 

DATA SWL,SWR/1,2/ 

Cc see ewer ee ewe Ee Ewe SB eee POS SH eB ee OR we OMe KH Bee eNO eK Bee Ce we wea wwwowen 
OPEN( UNI T<-20, FILE='WRPERF1.OUT’ ,STATUS='NEW’ )} 
WRITE(6,*)’ INPUT NUMBER OF TIME STEPS (2000):’ 
PEAD(S,*)K 
WRITE(20,*)' FILE WRPERF1.OUT FROM WRCOMPIA.FOR PROGRAM:’ 
WRITE(20,*) 

WRITE(20,*)!° $ OF TIME STEFS K= ‘',K 

WRITE(6,*)'INPUT CFL NUMBER: (RANGE 0.2 < CFLNUM <0.8)’ 

READ(5,*)CFLNUM 

WRITE(20,*%)}° CFLNUM= ',CFLNUM 7 
Cc Se nee Oo ok ee as 


WRITE(6,*)’ INPUT MAX NUMBER OF RIEMANN ITERATIONS (20):’ 
READ(S,*)QSTOP 
WEITE( 20,3021 )QSTOP 

3621 FORMAT(5X,'# OF RIEMANN ITER. (QSTOP)= ',13) 
WRITE(6,*)'INPUT GRID CELL WIDTH (0.01):° 
READ(5,*)DX 
WRITE(20,3022)DX 

3022 FORMAT(5X,'GRID CELL WIDTH DX= ',F4.2) 

C~----+-- INPUT TOTAL CONDITIONS INTO THE WAVE ROTOR----------~---- 
WRITE(6,*)' INPUT TOTAL PRESSURE IN (PTO) (psi):’ 
READ(5,*)PTO 

FTOTINePTO*144.%47.88026 


WRITE(6,*%)’ PTO =',PTO,'’psi’ 
WRITE(6,*)? PTOTIN@=’,PTOTIN, ‘Pa’ 
WRITE(20,*%)’ PTO =',PTO,’psi’ 
WRITE(20,*)’ PTOTIN=’,PTOTIN,'Pa’ 


WPITE(6,*)'INPUT TOTAL TEMP IN (TTOTIN) (R):° 

READ(S,*)TTOTIN 

WRITE(6,*)’ INPUT GAMMA IN:° 

READ(S,*)GAMMA 

G=GAMMA 

WRITE(6,*)* INPUT FUEL-AIR RATIO IN (FIN):° 

READ(5,*)FIN 
MGIN@|(1.4FIN)/(0.0345224+0.035648*FIN) 
RGIN@=1545.43/MGIN 


C-+-------- density in lbm/ft°3 
RTIN=©PTO*144.0/(RGIN#TTOTIN) 
C--------- density in kg/m°3 


RTOTIN@RTIN®0O. 453594(1. 70. 3048S FFT, 
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WRITE(6,*)° RTIN ~ ',RTIN,’ lbm/ft 3’ 
WRITE(6,*)’ RTOTIN = ',RTOTIN,’ kg/m™3’ 
WRITE(20,%)° RTIN = *,RTIN,’ lbm/£t°3" 
WRITE(20,*)° RTOTIN = ',RTOTIN,' kg/m*3' 


WRITE(6,*)*‘ INPUT TOTAL PRESSURE RATIO (TPR):’ 
READ(5,*) TPR 
WRITE(6,*)°INPUT EXIT MACH (ABSOLUTE):° 
READ(5,* )AMERIT 
C---~--~ COMPUTE EXIT STATIC PRESSURE 
PSEXIT1L=PTO*TPR/(14+(G-1.)/2.*AMEXIT*#*2,)**(G/(G-1.)) 
PSEXIT=PSEXIT1*144,0447,88026 


WRITE(6,*)° PSEXIT] = ',PSEXIT1,’psi' 
WRITE(6,*)° PSEXIT = ',PSEXIT,'Pa’ 
WRITE(20,¥*)" PSEXIT1 = ',PSEXIT1,’psi’ 
WRITE(20,")" PSEXIT ~ ',PSEXIT,'Pa’ 
C-----~- INPUT REF DATA-~--~----3~---~-3- 3-2 - eo woertee 


WRITE(6,4)'INPUT STATIC PRESSURE RATIO (SPR): i 
READ(S,*)SPR 


C----~-- COMPUTE PREF-~---~----~----~--~-+-~~--------~~---- 
PREFL*©PSEXIT1/SPR 
PREF=PREF1*144.0*47. 88026 
WRITE(6,*)}/ PREF1 = ',PREF1,’ (psi)’ 
WRITE(6,*)’ PREF «© ', PREF,’ (Pa)! 
WRITE( 20,4}! PREFL = ',PREFI,’ (pst)? 
WRITE( 22,4)! PREF «= °,PREF,’ (Pa)? 


WRITE(6,*) ‘INPUT GUESS FOR RREF (lbm/f£t*3):’ 
READ(5,*) RGUESS 
RREFeRGUESS*0, 45359*(1./0.304B6)**3. 


WRITE(6,*)/ RGUESS = ',RGUESS,' lbm/f£t*3! 
WRITE(6,*)! RREF «© ',RREF,' kg/m73! 
C~-----~- INPUT XREF IN METERS---~-----~---~~----------~-- 


WRITE(6,*) ‘INPUT XREF: (0.1800 m)’ 
READ(S,*)XREF 
C~------ CONVERT XREF TO FEET FOR CORRECT DIMENSIONS -~--~-- 
XREFL]@XREF/0O.3048 
WRITE(6,*)¢ xcef = ',XREF,' m' 
WRITE.:S,*)¢ xrefl = ',XREF1,’ £t’ 
WRITE( 20 )? xref = ',XREF,’ m’ 
WRITZ(?),*)? xrefi = ",XREFI,’ ft’ 


* 
* 


~~. 


GLU=GAMMA 
GR=GAMMA 
EPS#1.E-06 
Cc QSTOP=20 
N1]90 
JCOUNT=0 
KCOUNT=0 
UEXMAX=0. 
MASSINGO. 
MASOUT#-0. 
Cn ewe ee ee ee START COMPUTER TIME CLOCK-----~- ea s5 
RINTT<#0.0 
RCTIME«SECNDS(RINTT) 
OPEN(UNIT#15,FILE~'WRCOMP1.DAT’ , STATUS=' NEW’ } 
c Ox-0.01 
AREF@eSORT( PREF /RREF ) 
TIMRETe*XREF /AREF 
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RMU=SOQRT((G~-1.)/(G+1.)) 
X(1)=-0.5*Dx 
ZETA=WDP(1) 
X11@DX*(WDP(0)-0.5) 
DO 25 182,203 
X(1)9X(I-1)4+6.5*px 
25 CONTINUE 
bO 35 181,100 
XARRAY (1) =X(1*241) 
35 CONTINUE 
INITIAL DATA 
CALL INIT 
CALL INIT2L(PSEXIT, PREF, RREF) 
CALL INIT2R(PSEXIT, PREF, RREF) 
CALL INIT3L({PSINL,RINL) 
CALL INIT3R(PSINR, RINR) 
NONDIMENSIONALIZATION 
DO 30 I#1,203,2 
P(I)@-P(1)/PREF 
R(J)=R(1)/RREF 
U(1)=U(1)/AREF 
A(1)=#A(1) /AREF 
S(1)"ALOG(P(I)/R( 1) **G) 
30 CONTINUE 
CALL PLOT1(R) 
DO 40 J=1,K 
NeNe+1 
XI1=DX*(WDP(0)-0.5) 
QPRINT#N/50 
oT-100. 
vo 50 re1,203,2 
DTTeCFLNUM*DX/(2.*AMAX1(ABS(A(1)+U(1})),ABS(A(1)-U(1T))})) 
DT=-AMINIL(DTT,DT) 
50 CONTINUE 
TIOTAL=TTOTAL4DT 
TIME*TTOTAL* TIMREF 
XI=-XII 
vo 60 I=1,201,2 
PL=P(1) 
RL*=R(I) 
UL*u(1) 
PR=P(14+2) 
RReR(I+2) 
UR=U(I+4+2) 
XITEMPeXI 
IF(1.EQ.1) XI@ABS(XI) 
IF((1.£0.201).AND.(x1.GT.0.0)) XI=-XI 
CALL GLIMM(QSTOP,PSTAR,USTAR,ASTAR) 
XI1eXITEMP 
P(14+1)°PGLIM 
R(1+1)=RGLIM 
U(I4+1)*°UGLIM 
60 CONTINUE 
CALL GG(SWL,SWR,N,TTOTAL, TIME, UEXMAX, PTOTIN, PREF ,MASSIN, MASOUT) 
bo 70 1#1,201,2 
IF(XI.LT.0.) GOTO 80 
P(14+2)=P( 141) 
WRITE(6,*)'°P(", 142,’ de’, P(1+2) 
R(I42)eR( 141) 
WRITE(6,*)°R(’,142,')=',R(14+2) 
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U(14+2)<U(3T+1) 
WRITE(6,*)'U(',142,' Ju" ,U(I+2) 
A(14+2)=SQRT(G*P(14+2)/R(14+2)) 
WRITE(6,*)°A('’,3142,')=',A(I+2) 
S(1+2)=ALOG(P(1+2)/R(14+2}**G) 
WRITE(6,*}'S(',142,° )=’,S(I4+2) 
GOTO 70 
80 P(1)=P( 141) 
R(1I)"R(I4+1) 
U(1)=U(T+1) 
A(I)=SQRT(G*P(I)/R(I)) 
S(I)=ALOG(P(I)/R(1)**G) 
70 CONTINUE 


IF(SWL.EQ.1) CALL BCL1 

IF(SWL.EQ.2) CALL BCL2(PSEXIi,?REF) 
IF(SWL.EQ.3) CALL BCL3(PSINL,RINL, PREF, RREF) 
IF(SWL.EQ.4) CALL BCL4{(PTOTIN, RTOTIN, PREF, RREF) 
1F(SWL.EQ.5) CALL BCLS 

IF(SWR.EQ.1} CALL BCRI 

IF(SWR.EQ.2) CALL BCR2(PSEXIT, PREF) 
IF(SWR.EQ.3) CALL BCR3(PSINR,RINR, PREF, RREF) 
IF(SWR.EQ.4) CALL BCR4(PTOTIN, RTOTIN, PREF,RREF) 
IF({SWR.EQ.5) CALL BCRS5 

IF(SWL.EQ.4) MASSIN=MASSIN+R(3)*U(3)*DT 


IF(SWR.EQ.2) MASOUT=MASOUT+R(201)*U(201)*DT 
IF(N.EQ.(QPRINT*50)) CALL PLOT2(N,K) 
40 CONTINUE 
CALL ENDPL(0)} 
CALL DONEPL 


-~- = CALCULATE PRESSURES AND TEMPS FROM INPUT INFO~-~-~-~- 

PTE#PTO*TPR 
PE@PTE/((1.4(G-1.)/2. *AMEXIT**2.)**(G/(G-1.))) 
PREF=PE/SPR 
TREF@PREF*144,/( RGIN*RGUESS ) 
TE*TREF*(SPR**((G-1.)/G) ) 
TTE@TE*(1.4(G-1.)/2.*AMEXIT**2. } 

TSSs= CALC TOTAL TEMP RATIO FROM REF DENSITY GUESS~--------~ 
TTRieTPR**((G-1.)/G) 


WRITE(20,*) 
WRITE( 20,3024) 
WRITE(20,3025)SPR,TPR,TTRL, FIN, AMEXIT 
WRITE( 20,3026) PTO, TTOTIN,G,RTIN, RTOTIN 
WRITE( 20,3027) PTE, PE, PREF,TTE,TE 
WRITE(20,3028)TREF,RGUESS, RREF 
3024 FORMAT(2X,’WAVE ROTOR PERFORMANCE MAP ANALYSIS:‘) 
3025 FORMAT(5X,'SPR#’,F6.2,2X,'TPR =',F8.2,2X,’TTR =',F6.4,2%, 
>'FIN @°,F6.4,2X, “AMEXIT=’,F6.4) 
3026 FORMAT(5X,’PTO=',F6.2,2X,'TOTIN=',F8.2,2X,'GAMMA=',F6.4,2X, 
>*RTIN]',F6.4,2X,'RTOTIN’ ,F6.4) 
3027 FORMAT(5X,’PTE=’,F6.2,2X,’°PE =’,F8.2,2K,°PREF =’, F6.2,2X, 
>'TTE ©°,F7.1,1%,'TE =',F7.1) 
3028 FORMAT(4X,’TREF=',F7.1,1X, RGUESS=’,"7.5,2X,'RREF= ',F6.4) 
RCTIME*SECNDS{RINTT)-RCTIME 
WRITE(6,*) COMPUTER RUN TIME = ',RCTIME,’ secs,’ 
WRITE(15,*)’COMPUTER RUN TIME = ',RCTIME,’ secs.’ 
WRITE(20,*) 
WRITE(20,*)’COMPUTER RUN TIME = ',RCTIME,’ gece.’ 
CLOSE(UNIT#=15) 
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CLOSE(UNIT=20) 

STOP 

END 
SRRRARAEAEARRRARARAREERRRRARRRERRRRREREERRERERREARERERRERETHEHERED DEER 

SUBROUTINE GLIMM(QSTOP, PSTAR,USTAR,ASTAR) 

INTEGER Q,QSTOP 

REAL ML.MR,MLN,MRN 

COMMON/GLIMM1/PGLIM, RGLIM, UGLIM, PL,RL,UL,PR,RR,UR,AL,AR,GL,GR,EPS 

COMMON/GLIMM2/DT,DX,X1 

DATA Q,ML,MR/0,100.,100./ 

PSTAR=0.5*(PL+PR) 

COEFL=SQRT( PL*RL) 

COEFR=SQRT(PR*RR) 


ALPHA=1. 
c BEGIN GODUNOV ITERATION 
30 Q=-Q+1 
IF(PSTAR.LT.EPS) PSTAR#=EPS 
Cc COMPUTE NEXT ITERATION FOR ML AND MR 


MLN=COEFL* PHI ( PSTAR, PL) 
MRN=COEFR* PHI ( PSTAR, PR) 
DIFML=ABS(MLN-ML) 

DIFMR=ABS (MRN-MR) 
ML=MLN 
MR=MRN 
c COMPUTE NEW PSTAR 
PTIL=PSTAR 
PSTAR=(UL-UR+PL/ML4+PR/MR)/(1./ML4+1./MR) 
PSTAR=ALPHA* PSTAR+(1.-ALPHA) *PTIL 
IF(Q.LE.QSTOP) GOTO 10 
IF(ABS(PSTAR-PTIL).LT.EPS) GOTO 20 
c COMPUTE NEW ALPHA 
ALPHA=0.5*ALPHA 
Q=0 
IF((1.-ALPHA).LT.EPS) GOTO 20 
10 IF(DIFML.GE.EPS) GOTO 30 
IF(DIFMR.GE.EPS) GOTO 30 
END OF GODUNOV ITERATION; COMPUTE USTAR 
20 USTAR~(PL-PR+ML*UL4MR*UR) /(ML4MR) 
BEGIN SAMPLING FROCEDURE 
IF (XI.LT.USTAR*DT) GO TO 40 
RIGHT SIDE; SELECT CASE OF SHOCK OR EXPANSION 
IF (PSTAR.LT.PR) GO TO 50 
RIGHT WAVE IS A SHOCK WAVE 
WR=UR+MR/RR 
IF (XI.LT.WR*DT) GO TO 60 
c RIGHT OF RIGHT SHOCK CASE 
RGLIM*RR 
PGLIM=PR 
UGLIM#UR 
RETURN 
c LEFT OF RIGHT SHOCK CASE 
60 RGLIM=-MR/(USTAR-WR) 
PGLIM=PSTAR 
UGLIM=“USTAR 
RETURN 
Cc RIGHT WAVE IS A RAREFACTION WAVE 
50 CONST=PR/RR**GR 
RSTAR~(ESTAR/CONST) **(1./GR) 
ASTAR*SQRT(GR*PSTAR/RSTAR) 
AR=SQRT(GR*PR/RR) 
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IF (X1.GE.(USTAR+ASTAR)*DT) GO TO 70 
c LEFT OF RIGHT FAN CASE 
RGLIM@RSTAR 
UGLIM=USTAR 
PGLIM=PSTAR 
LTAG=3 
RETURN 
c SELECT RIGHT OF FAN OR IN FAN 
70 IF (XI.GE.(UR+AR)*DT) GO TO 80 
Cc IN RIGHT FAN CASE 
UGLIM=2./(GR4+1.)*(XI/DT-AR+(GR-1.)/2.*UR) 
RGLIM=(((AR4+(GR-1.)/2.*(UGLIM-UR) )**2.)/(GR*CONST) )**(1./(GR-1.)) 
PGLIM=CONST*RGLIM**GR 
RETURN 
c RIGHT OF RIGHT FAN CASE 
60 RGLIM=RR 
PGLIM=PR 
UGLIMaUR 
RETURN 
c LEFT SIDE; SELECT CASE OF SHOCK OR RAREFACTION 
40 IF (PSTAR.LT.FL) GO TO 90 
c LEFT WAVE IS A SHOCK WAVE 
WL=UL-ML/RL 
IF (XI.GE.WL*DT) GO TO 100 
Cc LEFT OF LEFT SHOCK CASE 
RGLIM=RL 
PGLIM]FL 
UGLIM=UL 
RETURN 
c RIGHT OF LEFT SHCCK CASE 
100 RGLIM=FL/(USTAR-WL) 
PGLIM=PSTAR 
UGLIM=USTAR 
RETURN 
c LEFT WAVE IS A RAREFACTION WAVE 
90 CONST#*FL/RL®* "GL 
RSTAR#=(PSTAR/CONST) **(1./GL) 
ASTAR@SQRT! Gu*PSTAR/RSTAR) 
AL=SQRT(GL*PL/RL) 
IF (XI.LT.¢(USTAR~ASTAR)*DT) GO TO 110 


c RIC 'T OF LEFT FAN CASE 
RGLIM=RSTAR 
PGLIM=PSTAR 
UGLIM=USTAR 
RETURN 
Cc SELECT LEFT OF FAN OR IN FAN CASE 
110 IF (XI,LT.(UL-AL)*DT) GO TO 120 
Cc IN LEFT FAN CASE 


UGLIM#=2./(Gu+i. )*(AL4+(GL-1.)/2.*UL4XI/0T) 
RGLIM=(((AL4+(GL-1.)/2.*(UL-UGLIM) )**2.)/(GL*CONST) )*#*(1./(GL-1.)) 
PGLIM=CONST*RGLIM* *GL 
RETURN 
Cc LEFT OF LEFT FAN CASE 
120 RGLIM#RL 
PGLIM=FL 
UGLIM*=UL 
RETURN 
END 
CORRARAEHHRAREARARE AER EERERERERAAARARERERRERARRERRERARAERAEEERERRREAD 


FUNCTION PHI(Y,2Z) 
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REAL RMU 

COMMON/FUNIL/G, PA, RA, UA, RB, RMU 

EPS=1.£-06 

PARAM=Y/Z 

IF (ABS(1.-PARAM).GE.EPS) Go TO 10 

PHI=SQRT(G) 

RETURN 

IF (PARAM.GE.1.) GO TO 20 

PHI@(G-1.)/2.*(1.—-PARAM) /(SQRT(G)*(1.-PARAM**((G-1.)/(2.*G6)))) 
RETURN 

PHI=SQRT((G+1.)/2.*PARAM+(G-1.)/2.) 

RETURN 

END 

FUNCTION PHI1(PB) 

REAL RMU 

COMMON/FUN1/G, PA, BA, UA, RB, RMU 

PHI l=(PB-PA) *SQRT((1.-RMU**2.)/(RA*(PB+RMU**2.*PA)) ) 
RETURN 

END 
RHRAERRHRARARARARARARAAARARARARARARRARARAARAAARARARRARARRHARRARARRHROERARHKEAN 
FUICTION PSI(PB) 

REAL RMU 

COMMON/FUNIL/G, PA, RA, UA, RB, RMU 
PSI#SQRT(1.-RMU**4.) /RMU**2./SQRT( RA) *PA**#(1./(2.*%G) )*( PB**((G~1.) 
>/(2.*G) )-PA**((G-1.)/(2.*G)})) 

RETURN 

ENS 

RR RARRAHRRARREKARRARARARARKRARRARARRRARAAHARARAAARHRAARKEHERAHRARARERRRE 
SUBROUTINE INIT1 

DIMENSION P(203),R(203),U(203),A(203),8(203),X(203) 
COMMON/FUN1/G, PA, RA, UA, RB, RMU 
COMMON/SUBS/P,R,U,A,S,X 

DO 10 r=1,9,2 

P(1)=-830600.00 

R(1I)=0.7132 

U(1)=644.4 

A(1)=SQRT(G*P(1I)/R(1I)) 

CONTINUE 

DO 20 1#11,203,2 

P(1)=101325.0 

R(I)=-1.22 

U(1)=0.0 

A(I)*SQRT(G*P(1)/R(I)) 

CONTINUE 

RETURN 

END 
RRAFARARAHRARARARHAKHARRARARAARRHRARARARAARAARARKRARRARRARKRAERKERRHREAL 
SUBROUTINE INIT2R( PSEXIT, PREF, RREF) 

DIMENSION P(203),R(203),U(203),A(203),5(203),X%(203) 
COMMON/FUN1L/G, PA, RA,UA,RB,RMU 
COMMON/SUBS/P,R,U,A,5,X 

DO 10 1re3,201,2 

P(I)=-PREF 

R(1)-RREF 

U(1)=-0.0 

A(I)-SQRT(P(I)*G/R(I)) 

CONTINUE 

P(1)=P(3) 

R(1)=R(3) 

U(1)=-u(3) 
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A(1)=SQRT(G*P(1)/R(1)) 
P(203)=PSEXIT 
R(203)=R( 201) 
PA=P(201) 
RA@#®R( 201) 
UA@U{ 261) 
PB=P(203) 
RB=R(203) 
IF(PA.GT.PB) GO TO 20 
U(203)=sUA-PHI1( PB) 
GO TO 30 
20 U(203)—UA-PSI(PB) 
30 A(203)=SQRT(G*P(203)/R(203)) 
RETURN 
END 
CAAARRHKHEERHRAARRRAREAHARAAEARAHRRARRARARHARARHRRARARRARARARRARE 
SUBROUTINE INIT2L(PSEXIT) 
DIMENSION P(203),R(203),U(203),A(203),8(203),X%(203) 
COMMON/FUNL/G, PA, RA, UA, RB, RMU 
COMMON/SUBS/P,R,U,A,S,X 
DO 10 183,201,2 
P(1)=285080.0 
R(1)=0.697 
U(1)=0.0 
A(1)-SQRT(G*P(I)/R(I)) 
10 CONTINUE 
RETURN 
END 
CARRRARKRHRAARHKANHKARARERHRHARRAARRAREHRARAARRAARREAHARARRERAARA 
SUBROUTINE INIT3L(PSINL,RINL) 
DINENSIC™ P(203),R(203),0(203),A(203),8(203),X%(203) 
COMMON/FUNL/G, PA, RA, UA, RB, RMU 
COMMON/SUBS/P,R,U,A,5,X 
pO 10 1=3,201,2 
P(r)=2390000.0 
R(1)=9.767 
u(1I)=0.0 
A(T) @SQni(G*P(I)/R(I}} 
10 CONTINUE 
Pll)sFETVL 
R(1L)*RINL 
PA=P( 3) 
RA#R(3) 
UA#U(3) 
PB=P(1) 
U(1)=sUA4+PHI1( PB) 
A(1)-SQRT(G*P(1)/R(1)) 
P( 2932 =P(201) 
R(203)=R( 201) 
U(203)8-u( 201) 
A(203)@SQRT(G*P(203)/R(203)) 
RETURN 
END 
CHARARHHAPHEAHARRARERARAHHARARERAHRARRERARARARKARHARARERARAEHRERS 
SUBROUTINE INIT3IR(PSINR, RINR) 
DIMENSION P(203),R(203),U(203),A(203),S8(203),X(203) 
COMMON /SUBS/P,R,U,A,S,X 
COMMON/FUN] /G, PA, RA,UA, RB, RMU 
DO 10 1=3,201,2 
P(1)#2421667.5 
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R(I)=9.787 
U(T)=0.0 
A(I)=SQRT(G*P(I)/R(I)) . 
10 CONTINUE 
P(203)=PSINR 
PA=P(201) 
RA=R( 201) 
UAeU( 201) : 
PB=P(203) 
U(203)=<UA~PHI1(PB) 
R(203)=RINR 
A(203)=SORT(G*P(203)/R(203)) 
P(1)=P(3) 
R(1)=R(3) 
U(1)=-u(3) 
A(1)=SQRT(G*P(1)/R(1)) 
RETURN 
END 
CARAAAAERHRAERARHARERARARRARRARRARARARRRARARHARERARAKERERAREREHEREH 
SUBROUTINE BCL1 
DIMENSION P(203),R(203),U(203),A(203),6(203),X(203) 
COMMON/SUBS/P,R,U,A,S,X 
COMMON/FUNI/G, PA, RA, UA, RB, RMU 
P(1)=P(3) 
R(1)=R(3) 
U(1)=#-U(3) 
A(1)=SQRT(G*P(1)/R(1)) 
RETURN 
END 
Creraank 
SUBROUTINE BCR3(PSINR,RINR, PREF, RREF) 
DIMENSION P(203),R(203),U(203),A(203),8(203),X(203) 
COMMON /SUBS/P,R,U,A,S,X 
COMMON/FUN1/G, PA, RA,UA, RB, RMU 
F(203)=PSINR/PREF 
R(203)@RINR/RREF 
PA=P(201) 
RA#R( 201) 
UA=U(201) 
PB=P( 203) 
U(203)=UA-PHI1(PB) 
A({203)=SQRT(G*P(203)/R(203)) 
RETURN 
END 
Crean 
SUBROUTINE BCL3(PSINL,RINL, PREF,RREF) 
DIMENSION P(203),R(203),U(203),A(203),8(203),X%(203) 
COMMON/SUBS/P,R,U,A,S,X 
COMMON/FUNI1/G, PA, RA, UA, RB, RMU 
P(1)=PSINL/PREF 
R(1)“RINL/RREF 
PA=P(3) 
RA#R( 3) 
UAeU(3) 
PB=P(1) 
U(1)sUA+PHI1(PB) 
A(1)@SQRT(G*P(1)/R(1)) 
RETURN 


END 
Cranes 
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Cesane 


Crekaan 


10 
20 


SUBROUTINE BCR1 

DIMENSION P(203),R(203),U(203),A(203),5(203),%(203) 
COMMON/SUBS/P,R,U,A,S,% 
COMMON,’FUN1/G, PA, RA, UA, RB, RMU 

P(263)=P( 201) 

R(203)"R(201) 

U(203)=-0(201) 

A(203)=SORT(G*P(203)/R(203)) 

RETURN 

END 


SUBROUTINE BCL2(PSEXIT, PREF) 

DIMENSION P(203),R(203),U(203),A(203),8(203),%(203) 
COMMON/SUBS/P,R,U,A,8,X% 
COMMON/FUNI/G, PA, RA, UA, BB, RMU 

P(1)<PSEXIT/PREF 

R(1)=R(3) 

ul )jeu(3) 

ACL )=SORT(G*P(L)/R(L)} 

RETURN 

END 


SUBROUTINE BCR2(PSEXIT, PREF) 
DIMENSION P(203),R(203),U(203)_A(203),5(203) ,x(203) 
COMMON/SUBS/P,R,U,A,S,% 
COMMON/FUNI/G, PA, RA, UA, RB, RMU 
P(203)=PSEXIT/PREF 
R(203)=R(201) 

PA~P(201) 

RA~R/201) 

VA=U( 201) 

PBeP( 203) 

RBeR( 203) 

If (PA.GT.FB) GO TO 10 
U(203)=#UA-PHI1( PB) 

GO TO 20 

U(203)=UA-PSI(PB) 
A(203)=SQRT(G*P(203)/R(203)) 
RETURN 

END 


CHAAHAARARHRAEAARREHEEHAREHERARRARHRERRRAHRERARERRHRHARARAE 


20 


10 


FUNCTION WDP(IT) 

DIMENSION WNORM(12),IDIGT(12) 
COMMON/SAMPLE/WNORM, IDIGT 
If (11.£Q.0) Go To 10 

Lle=2 

L2=1 

bo 20 JJe1,12 

IDIGT(JI)}=0 
WNORM(JJ)=1./FLOAT(LI**JIJ) 
CONTINUE 

WDPeD. 

RETURN 

DO 40 JJe1,12 

L1-2 

L2=1 

RID“IDIGT(JII) 

RIN=MOD({ (KJ0+1),L1) 
IDIGT(JJ) KIN 

le (RIO.LT.KIN} GO TO 50 
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40 CONTINUE 


50 SuUM=0. 


bo 60 


JJdu1,12 


KNEW=MOD(IDIGT(JJ}*L2,01) 

SUM=SUM+ FLOAT (KNew) *WNORM(JJ) 
60 CONTINUE 

WOP=SUM 


RETURN 


END 


CHARARARRAAHRHEHRARAARAEAREREARKARRARAARERRERRAERARRARHARRARARARRKARRERARS 


>00) 


DATA 
DATA 
DATA 
DATA 
DATA 
DATA 


CALL 
CALL 
CALL 
CALL 
CALL 


IF(I 
Ie(? 
IF(I 


IF(I 
IF(I 
IF(I 
CALL 


END 


DATA 
DATA 
DATA 
DATA 


CALL 
CALL 


ANANANAANANAANANANANAAANAANAANAANAAANANANANNANANNANANANANANANANNNNANNANNANANANA 


IF(I. 


bo 10 


SUBROUTINE PLOT2(N,K} 
DIMENSION XORG(4),YORG(4),YMAX(4),YMIN(4),RNT(4),1YNAM(10) 
DIMENSION PARRAY(100),RARRAY({100),UARRAY(100),SARRAY(100),XARRAY(1 


DIMENSION P(203),R(203),U(203),A(203),5(203) ,X(203) 
COMMON/SUBS/P,R,U,A,S,X% 
COMMON XARRAY 


XORG/O.5,4.75,0.5,4.75/ 

YORG/0.5,0.5,4.75,4.75/ 

¥MAX/2.50,2.5,2.5,2.5/ 

YMIN/0.0,0.0,-1.0,-2.0/ 

KNT/1,4,5,9/ 

IYNAM/'PRES','SURE’,’S ', DENS’, ‘ITY$',’VELO’,'CITY',’$ : 


>, "ENTR’, 'OPY$'/ 
DO 200 1,100 
PARRAY(1)#P(1*241) 
RARRAY(1)=R(1*2+1) 
UARRAY(1)=U(1*2+1) 
SARRAY(I)=S(1*2+1) 
200 CONTINUE 
DO 300 11,4 


PHYSOR( XORG(1I), YORG(1I)) 

AREA2D(3.5,3.5) 

XNAME('X',1) 

YNAME(IYNAM(KNT(1I)),100) 
GRAF(0.,’SCALE’,12.0,YMIN(I),'SCALE',YMAX(I)) 
£Q.1) CALL SETCLR(’ YELLOW’ ) 


-EQ.2) CALL SETCLR( ‘CYAN’ ) 
-EQ.3) CALL SETCLR(’RED’) 
-EQ.4) CALL SETCLR('MAGENTA’ ) 
IF(N. 
c IF(N. 
IF(I. 


£Q.K) CALL SETCLR('WHITE’ ) 
EQ.2398) CALL SETCLR(’BLUE’ ) 
EQ.1) CALL CURVE (XARRAY, PARRAY,100,0) 


-EQ.2) CALL CURVE (XARRAY,RARRAY,100,0) 
-EQ.3) CALL CURVE (XARRAY,UARRAY,100,0) 
.EQ.4) CALL CURVE (XARRAY,SARRAY,100,0) 


ENDGR(0) 


300 CONTINUE 
RETURN 


PREREHRRERRAEEAAER EER EHREREREARERERERRRA ARR RRARRARREREAHED 
SUBROUTINE PLOT1(K) 
DIMENSION XORG(4),YORG(4),¥YMAX(4),YMIN(4) 


XORG/0.5,4.75,0.5,4.75/ 

YORG/O.5,0.5,4.75,4.75/ 

Y¥MAX/2.50,2.5,2.5,2.5/ 

¥MIN/O.0,0.0,-1.0,-2.0/ 
io1,4 

PHYSO®?(XORG(I), YORG(I)) 

AREA2D(3.5,3.5) 
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CALL FRAME 
CALL GRAF(0.,’SCALE’,1.0,YMIN(I),‘SCALE’ , YMAX(1)) 
CALL ENDGR(0) 
10 CONTINUE 
CALL PHYSOR(6.5,0.5) 
CALL AREA2D(2.25,7.75) 
CALL FRAME 
CALL GRAF(0.,'’SCALE’,1.,0,'SCALE’ ,K) 
CALL ENDGR(0) 
RETURN 
END 
RRARKARAMAARAARRARRARARARKRARAARAKRAAAARAAKRKHEARREARHRAAAREARARKAAARER 
SUBROUTINE BCL4 (PTOTIN,RTOTIN, PREF, RREF) 
INTEGER QOUT 
DIMENSION P(203),R(203),U(203),A(203),8(203),X%(203) 
DIMENSION QR(3) 
DIMENSION XARRAY(100) 
COMMON/SUBS/P,R,U,A,S,X 
CONNOI/FUNL,/G,PA,RA,UA,RB, RMU 
COMMON XARRAY,N1 
NileN141 
QR(3)"2./(G-1.)#*A(3)-U(3) 
QOUT=N1/5 
PTOT#FIOTIN/PREF 
RTOT@RTOTIN/RREF 
ATOT=SQRT(G*FTOT/RTOT) 
STOT=“ALOG(PTOT/RTOT**G) 
u(1seu(3) 
A(1)=SQRT(ATOT**2,-(G-1.)/2.*ABS(U(1))**2.) 
OR(1)~2.,/(C-1.)*A(1L)-u(1) 
AMACH®U(1)/A(1) 
IF(AMACH.LT.0.0) GO TO 60 
P(1)=PTOT/(1.+(G-1.)/2.*AMACH**2.)**(G67(G-1.)) 
R(1)=RTOT/(1.4(G-1.)/2.*AMACH*#2.)**(1./(G-1.)) 
§(1)-ALOG(P(1)/R(1)**G) 
GO TO 50 
60 P(l)=PTOI/(1.+(G-1.)/2.*ABS(AMACH)**2.)**(G/(G-1.)) 
R(1)@-RTOT/(1.4+(G-1.)/2.*ABS(AMACH)**2.)**(1./(G-1.)) 
S(1L)*ALOG(«(1)/R(1)**G) 
50 TF(N1.€Q.5*QOUT) WRITE(10,1) A(1),A(3),U(1),0(3),N1 
1 FORMAT(2™%,4F9.6,5X,15) 
7 ‘TURN 
ED 
CH anah 
SUBROUTINE BCR4 (PTOTIN,RTOTIN, PREF,RREF) 
INTEGER QOUT 
DIMENSION P(203),R(203),U(203),A(203),8(203),%(203) 
DIMENSION QL(203) 
DIMENSION XARRAY(190) 
COMMON /SUSS/P,R,U,A,S,X 
COMMON/FUN1/G, PA, RA, UA, RB, RMU 
COMMON XARRAY,N1 
Nl=Nl41 
QL(291)"2./(G-1.)*A(201)4U(201) 
QOUTHNi/25 
PTOT=PTOTIN/PPEF 
RTOT#RTIOTIN/RREF 
ATOT=SORT(G*PTOT/RTOT) 
STOT#ALOG( FTOT/RTOT* #G) 
U(203)©U(201) 


AANAAAAANAANNA 





60 


Creare 


Crete 


Crake 


10 


20 


30 


A(203)*SQRT(ATOT**2.-(G-1.)/2.*ABS(U( 203) )**2.) 
QLU(203)—2./(G-1. J) *A(205)4U( 203} 

AMACH=U( 203) /A( 203) 

IF(AMACH.LT.0.0) GO TO 60 
P(203)-PTOT/(1.+(G-1.)/2.*AMACH* #2.)#*(G/(G-1.)) 
R(203)-RTOT/(1.4(G-1.)/2.*ANACH**2.)#*(2./(G-1.)) 
$(203)-ALOG(P(203)/R(203)**G) 

GO TO 50 

P(203)@PTOT/(1.4(G-1.)/2.*ABS(AMACH) #*2.)**(G/(G-1.)) 
R(203)RTOT/(1.4(G-1.)/2.*ABS(AMACH) #*2.)**#(2./(G-1.)) 
§(203)-ALOG(P(203)/R(203)**G) 

LECNL.EQ.25*Q00UT) WRITE(10,1) P(201),0(201),N1L,N 
FORMAT(2X,2F13.6,5X,215) 

RETURN 

END 
* 

SUBROUTINE BCL5 

DIMENSION P(203),R(2063),0U(203),A(203),8(203),X%(203) 
COMMON,/SUBS/P,R,U,A,5,X% 

p(1)=P( 3} 

R(1)~R(3) 

u(1)-U(3) 

A(1)=A(3) 

RETURN 

END 
* 

SUBROUTINE BCRS 

DIMENSION P(203),R(203),0(203),A(203),8(203),X(203) 
COMMON/SUBS/P,R,U,A,8,X 

P(203)=sP{201) 

R(203)=R(201) 

U(203)=U(201) 

A(203)=8A(201) 

RETURN 

END 
KERARRAARARHERARRAAAAARAHRRAAERNRERRAERRAAKRHAARARARAAHAHHHERKRAKARARA RA RAH 
SUBROUTINE CC(SWL,SWR,N,TTOTAL, TIME, UEXMAX, PTOTIN, PREF ,MASSIN,MASO 
*uUT) 

INTEGER SWL,SWR 

R. \L MASSIN, MASOUT 

DIMENSION P{(203),R(203),U(203),A(203),8(203),X%(203) 
COMMON/SUBS/P,F,U,A,8,X% 
IF((SWL.EQ.1).AND.(SWR.EQ.2)) GO TO 19 
IF((SWL.EQ.4).AND.(SWR.EQ.2)) GO TO 30 
IF((SWL.EQ.4).AND.(SWR.EQ.1)}) GO TO 50 
IF((SWL,EQ.1).AND.(SWR.EQ.1)) RETURN 

PWALLeP(2) 

IF(PWALL.LE.(PTOTIN/PREF}) GO TO 20 

RETURN 

SWLe4 

WRITE(6,74) 

WRITE(6,75) N,TTOTAL, TIME 

WRITE(15,74) 

WRITE(15,75) N,TTOTAL, TIME 

WRITE(20,*%) 
WRITE(20,*)’TIMING AND MASS FLOW BALANCE:' 

WRITE(20,74) 

WRITE(20,75) N,TTOTAL, TIME 

RETURN 

VUEXIT=U( 202) 
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40 


50 


60 


73 
74 
75 
78 
76 
77 


IF (UEXMAX.LT.UEXIT) UEXMAX*UEXIT 
IF(UEXIT.LT.UEXMAK/2.) GO TO 40 
RETURN 
SWR-1 
WRITE(6,76) 
WRITE(6,75) N,TTOTAL, TIME 
WRITE(15,76) 
WRITE(15,75) N,TTOTAL, TIME 
WRITE(20,76) 
WRITE(20,75) N,TTOTAL, TIME 
RETURN 
PISHOR=P(2) 
IF(N.EQ.1590) GO TO 60 
IF(PISHOK.GT.PTOTIN/PREF) GO TO 60 
RETURN 
SwLeal 
WRITE(6,77) 
WRITE(15,77) 
WRITE( 20,77) 
WRITE(6,75) N,TTOTAL, TIME 
WRITE(15,75) N,TTOTAL, TIME 
WRITE( 20,75) N,TTOTAL, TIME 
WRITE(6,78) P(199),P(201),P(203) 
WRITE(15,78) B(199},P(201),P(203) 
WRITE(20,78) P(199),P(201),P(203) 
WRITE(6,73) MASSIN,MASOUT,SWL,SWR 
WRITE(15,73) MASSIN,MASOUT,SWL,SWR 
WRITE(20,73) MASSIN,MASOUT, SWL,SWR 


FORMAT(SX,'Mine’,F9.6,2xX, ‘Mout=’,F9.6,5X,'SWL='’,14,2X,'SWRe’ 


FORMAT(SX,’INLET PORT OFENS AT:"') 


FORMAT(5X,'N=’,14,5X,*TTOTAL@',F14.7,2X,'TIME@’,F14.7) 


14) 


FORMAT(lX,'P(199})=°,F8.6,2K,'P(201)=',P6.6,2X,°P(203)=',F8.6) 


FORMAT(5X,'EXHAUST PORT CLOSES AT:’) 
FORMAT(5X,’ INLET PORT CLOSES AT:’) 
RETURN 

END 
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APPENDIX C 


COMPUTER SYSTEM AIDS 


TRANSFERRING FILES 
1. FROM MAINFRAME TO AERO PC LAB 
The following procedures use SIMPC to download 
programs from the NPS mainframe to the Aero Department's PC 
Lab. The modem and SIMPC hooked up to the IBM PC-XT-2 in the 
Aero PC lab. The following boldfaced commands outline the 
procedures. 


C\ <dir> >cd\SIMPC......... change dir to SIMPC dir 
C\SIMPC>SIMPC ACCESS.SIM...execute SIMPC macro ACCESS.SIM 


LOGON: 9812P ...-logon to mainframe (user number) 
PASSWORD: --..enter your mainframe password 





<enter> +...@xecutes your PROFILE.EXEC 


FORSIMPC GET ENGINE FORTRAN Al..file to copy from 
mainframe 


A: ENGINE.FOR 


FORSIMPC PRINT ENGINE FORTRAN Al..optional if you wanted 
to send your file to the current printer. 


LOGOFF -.-logoff the mainframe 
<ALT> <G> -.e-hangs up the modem 
<ALT> <Q> -.-quits SIMPC and returns to DOS 
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The procedure above uses the macro ACCESS.SIM which 


was written specifically to link between the Aero Pc lab nd 
the NPS mainframe. 

2. FROM AERO P.C. LAB TO VAX SYSTEM 

The following procedure uses PROCOMM, an executable 

modem link, and ALT 0, a modem macro written specifically to 
link the Aero Pc lab to the Aero VAX lab. KERMIT is the 
actual file transferring program. 

C\ <dir> >cd\COM...change to COM directory. 


FROCGIM .executcble program 
<ALT> O .macro modem link 


ATDT2953 

CONNECT 

USERNAME: GUEST -gouge username 

PASSWORD: -password will not appear 


Are you logged on from a PC? Y 
Do you want to connect to WASP? Y 


Username: -typical VAX logon 
Password: . 
C\\>KERMIT ..executes file transfer program 


to send a file: 

Receive *RCOMP1.FOR...Program you wish to call 

F3: WRL ».FOR’ ... file you wish to send but this file 
must be in the PROCOMM dir on the C (hard) drive! 


LOGOFF -.--logoff the VAX 
<ALT> x ...exits PROCOMM 


If you have any difficulties F10 is the help menu! 
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B. USING "SIGHT" UTILITY ON THE VAX 

SIGHT is a utility drawing program available on the Aero 
Department VAX computer. It was used to create pop-up menus 
for the configurations for the ENGINE program. To run SIGHT, 
at the DCL prompt enter 

SSIGHT 
A pop-up menu provides simple instructions. Save your drawing 
as a <fn>.UIS file and then use the VAX "render" command to 
create your output file as a <fn>.REN file. An example for 
sending creating and sending an output file to the laser 
(1n03) printer is: 

SRENDER <fn>.UIS/DEVICE TYPE=1n03 

$PRINT/PASSALL/QUE=1n03 fn.REN 
SIGHT can be used to generate a picture from inside a program. 
This can be done by inserting the following statement in a 
(FORTRAN) program. 


CALL LIBSSPAWN ('RENDER <fn>.UIS") 


This was done for the ENGINE program. 





C. GRAFkKit PROCEDURES 

PLOT3.FOR is a GRAFkit plotting program which was 
incorporated in the ENGINE1A.FOR program. PLOT3.FOR will be 
available as an engine performance plotting program on the VAX 
network. It can accept up to 50 separate data inputs for a 
selected X-axis, a Yl-axis for specific thrust, and a Y2-axis 
for specific fuel consumption. The instructions for running 
PLOT3.FOR are included in the comments section of the program 


listing. 


PLOT3.FOR program listing 


THIS PLOTTING PROGRAM USES GRAFKIT SOFTWARE 
ENSURE THE FOLLOWING SEQU2NCE IS ENTERED PRIOR TO 
EXECUTION OF PLOT3: 
Isgon the VAX 
GRSETUP 
$urs 
$RUN PLOT3 
If YOU WANT YOUR RESULTS SENT TO EITHER THE LASER OR 
HIGH SPEED PRINTERS: 
logon the VAX 
SGRSETUP 
$1n03 
$DEFINE GRK_OUT PLOT.DAT ("PLOT.DAT" IS YOUR) 
$1n03s8_ep (OUTPUT fn.ft ) 
SRUN PLOT3 
«eto send output to the selected printer: 
SFRINI/PASSALL/QUE=LNO3 PLOT.DAT 


c 
Cc 
Cc 
Cc 
Cc 
c 
Cc 
Cc 
c 
Cc 
c 
c 
Cc 
c 
c 
Cc 
Cc 


PROGRAM PLOT3 

REAL X1(50),¥1(50),¥2(50) 

CHARACTER X1LABEL*10, YILABEL* 33, Y2LABEL* 33 
CHARACTER FIG1*7,NUMB1*3,TOPLBL*10 


Camm ewe were eee nes ane sree nals Se senses sennesenaeeneeneene 


SPECIFY X,Y LEFT AND Y RIGHT LABELS 


YLULAREL@’A: Specific Thrust Ibf£/(lbm/e}" 
Y2ZLABEL~’B: SFC {lbm fuelyhr)/ibf thrust’ 
CHOOSE THE X PARAMETER TO PLOT 
WRITE(6,*%)’ CHOOSE YOUR X AXIS LABEL({i.e. 1):’ 
WRITE(6,*)" 1. BYPASS RATIO’ 
WRITE(6,*)’* 2. FAN PRESSURE RATIO’ 
WRITE(6,*)° 3. COMPRESSOR RATIO’ 
WRITE(6,*%)° 4. LPT EXIT MACH NUMBER’ 
WRITE(6,*)° 5. OTHER (MAX 10 CHARACTERS)’ 
READ(5,*)IANSS 
IF (IANS4.£Q.1) THEN 
XILABEL@' BETA ' 
ELSEIF (1ANS4.£Q.2) THEN 
X1LABEL=’ PRLC ’ 
ELSEIF (IANS4.£Q.3) THEN 
XILABEL=* PRHC * 
ELSEIF (1ANS4.£Q.4) TREN 
XILABEL@='AMACHB’ 
ELSEIF (TANS4.£Q.5) THEN 
WRiIE(6,4)° INeul AX AAIS LAEE:! 
READ(S,°(A10)')XILABEL 
ENDIF 
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C--.---CHOOSE THE NUMBER OF DATA PTS TO PLOT MAK=50 


100¢ 


Cc 


Cewek 26 ee 8 E688 8S 8 BES SO EES SSE EERE SESeRSSes Ses ARSENE 








WRITE(6,*)' INPUT NUMBER OF DATA POINTS TO PLOT (MAX=50)' 
READ(S,*) ICP) 
INPUT X AXIS VALUES: 


WRITE(6,*)' INPUT VALUES FOR ',X1LABEL 
DO rel, Icel 
WRITE(6,1000)°X(",3,') @! 
READ(5,*)X1(T) 
END DO 


WRITE(6,*)° INPUT? 
WRITE(6,*)?’ 0 - AUTOMATICALLY DETERMINE X1MIN AND X1MAX’ 
WRITE(6,%)° 1 - INPUT X1MIN AND X1MAX’ 
READ(5,*) IANS5S 
If (IANS5.EQ.1) THEN 
WRITE(6,*)* INPUT XIMIN =3’ 
PEAD(5,*)X1MIN 
WRITE(6,¢)* INPUT XIMAX o:’ 
READ(S,*)X1MAX 
AUTOMATICALLY DETERMINE THE 
ELSE 
XIMIN@=X1(1) 
XIMAX=K1(1) 
DO Jel, ICPl 
XITEMP=x1(9) 
X2TEMPHX1(9) 
IF {X1MIN.GT.XITEMP) THEN 
XIMIN@XLTEMP 
ENDIF 
IF (XIMAX.LT.X2TEMP) THEN 
KIMAReX2TEMP 
ENDIF 
£ND DO 
ENDIF 
INPUT VALUES FOR YL AND ¥2 ARRAYS 
WRITE(6,*)' INPUT VALUES FOR ‘',YILABEL 
DO Jeli, ICcPl 
WRITE(6,1000)°¥1(',9,') = 
READ(5,*)¥2(9) 
END DO 
WRITE(6,*)' INPUT VALUES FOR ',Y2LABEL 
DO K=1,1CP1 
WRITE(6,1000)’¥2(",K,’) =? 
READ(5,*)} Y2(K) 
END 0O 
FORMAT(1X,A3,12,A3) 


MIN AND MAX K VALUE 


INPUT THE NUMBER OF THY FIGURE TO PLOT (FIG 1)--- 
WRITE(6,*)'INPUT TRE FIGURE NUMBER (fe 1 OR 2...)' 
READ(5,°(A9)’ )NUMBL 
WRITE(6,*)"NUMB1e’ ,NUMB1] 
FIGL='FIGURE ° 
TOPLBLe@FIG1//NUMB} 
WRITE(6,°(1K,A11)’)* TOP LABEL =’, TOPLBL 


CONSTRUCT THE PLOT 


WRITE(6,*)’STANDBY FOR PLOT’ 
OPEN GKS, OPEN AND ACTIVATE WORKSTATION. 
CALL OPNGKS(1,2,1) 


SET UP 
CALL 
CALL 
CALL 


THE COLOR TABLE 

Gscr(1,0,0.00,0.00,0.50) 
GScR(1,1,1.00,1.00,0.00) 
GsScr(1,2,0.00,1.00,0.50) 


CALL GSCR(1,3,0.09,1.00,1.00) 

CALL GSCR(1,4,1.00,0.00,0.00) 

CALL GSCR(1,5,1.00,0.00,1.00}) 

CALL GSCR(1,6,1.00,0.00,1.00) 

CALL G8CR(1,7,1.00,1.00,1.00) 
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CALL GSCR(1,8,0.00,0.00,0.00) 


c 
C----- SET VIEWPOINT AND DEFINE THE ad USER'S DATA TO BE USED 

CALL AGSETF ('GRID/LEFT.’, 1 

CALL AGSETF ('GRID/RIGHT.’, .85) 

CALL AGSETF (‘GRiD/TOP.’, .85} 

CALL AGSETF (’GRID/BOTTOM.’, .15) 
C----- SET UP FOR FIRST PLOT 
C----- DEFINE THE TOP, BOTTOM, LEFT AND RIGHT LABELS 
C----- PUT BLANKS IN FIRST CURVE LABELS 

CALL AGSETC (’LABEL/NAME.’, ‘T') 

CALL AGSETI ('LINE/NUMBER.*, 100) 

CALL AGSETC (’LINE/TEXT.', ee | 
c 

CALL AGSETC (‘LABEL/NAME.’, ‘B’} 

CALL AGSETI (’LINE/NUMBER.',-100) 

CALL AGSETC ('LINE/TEXT.’, ef) 
c 

CALL AGSETC ({'LABEL/NAME.’, ‘L’) 

CALL AGSETI (’LINE/NUMBER.’, 100) 

CALL AGSETC ('LINE/TEXT.', rot) 

CALL AGSETC (’LABEL/NAME.', ‘'R') 

CALL AGSETI ('LINE/NUMBER.’,-100)} 

CALL AGSETC (’LINE/TEXT.', ery 
C----~ TURN OFF RIGHT AXIS 

CALL AGSETF (’AXIS/RIGHT/CONTROL.’ ,0.0) 
C----- SET YMIN AND YMAX FOR FIRST PLOT 


CALL AGSETF (’K/MINIMUM.%, X1MIN) 
CALL AGSETF ("X/MAXIMUM.‘', KIMAX) 
CALL AGSETF ('Y/MINIMUM.’, 0.0) 
CALL AGSETY ('Y/MAXIMUM.’, 120.0) 


CALL AGSTUP (X1,1,1,ICP1,1,¥1,1,1,1CP1,1) 

CALL AGBACK 

CALL GSPLCI (7) 

CALL DASHDC ('SSSSSSSSSSS' ITT hr SSssssssss’,3,2) 
C------ PLOT CURVE 1: xl VS YI 

CALL AGCURV (X1,1,¥1,1,ICP1,-1) 


Coaessesveseanesussesenes 
Cases SET UP FOR SECOND PLOT 
Caen DEFINE THE TOP, BOTTOM, LEFT AND RIGHT LABELS 
Connn- THE LARGER LINE NUMBER WRITES ABOVE PREVIOUS LABELS 
CALL AGSETC ('LABEL/NAME.'’, Tv’) 
CALL AGSETI ('LINE/NUMBER.’, 100) 
CALL AGSETC (’LINE/TEXT.'’,TOPLBL) 


CALL AGSETC (‘LABEL/NAME.’, ‘'B’) 
CALL AGSETI ('LINE/NUMBER.*,-110) 
CALL AGSETC (’LINE/TEXT.’, KILABEL) 


CALL AGSETC (‘LABEL/NAME.’,- "L’) 
CALL AGSETI (’LINE/NUMBER.’, 110) 
CALL AGSETC (*LINE/TEXT.’ , YILABEL) 


CALL AGSETC ("LABEL/NAME.', = 'R') 
CALL AGSETI ('LINE/NUMBER.’, -100) 
CALL AGSETC (’LINE/TEXT.', Y2LABEL) 

c TURN OFF THE TOP, BOTTOM AND LEFT AXES 
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eee ee eee =- oon” a= 


CALL AGSETE (SARIS/TOP/CONTROL.' 10.0) 
CALL AGSETF (AKIS/BOTTOM/CONTROL. *,0-0) 
CALL AGSETFE (tAKIS/LEFT/CONTROL.',0-0) 


C-~--" TURN ON THE RIGHT AXI 


3S 
CAL. AGSETY ('RR1S/RIGHT/CONTROL. ', 3-0) 


C-~---~ CRANGE THE YMIN AND YMAX 


CALL AGSETE (7 Y/MINIMUM. ', 0.0) 
CALL AGSETF ('Y/MAXIMUN.” y 1.5) 
C----- TURN ON THE RIGHT AXIS NUNERICS 
CALL AGSETE ( RIGRT/NUMERIC/TYPE.’ ,1-¥6) 


c 

CRLL AGSTUP (X1,4,2,1C 

CALL AGBACK 

CALL GSPLCI (4) 

CALL DASHOC U'gssEcs'’ 
C--0077 PLOT CURVE 21 Xl vs ¥2 

CALL AGCURV (R1,2,¥2,4 
C--~~-~ CALL FRAME TO ADVANCE TH 
Coorre CLOSE GXS, CLOSE AND DEA 


CALL FRAME 
CALL CLSGRS 


Caseewr canes ceussavesnaee 


£ND 


94 





PLeL,Y¥2,helsICPL-2) 


e+ 8,1) 


,ICPl,-2) 
& FRAME 
CTIVATE WORKSTATION 





10. 


11. 


le. 
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